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Executive  Summary 


Digital  wetland  mapping  provides  important  infor- 
mation on  the  t\'pe,  location,  and  extent  of  wetlands 
within  a  given  region.  Comparing  historic  map- 
ping with  updated  mapping  provides  a  unique  op- 
portunity to  examine  potential  changes  in  wetland 
density  and  distribution  due  to  both  natural  and 
anthropogenic  causes.  In  addition  to  documenting 
changes  in  wetland  area,  comparing  spatial  datasets 
allows  us  to  track  change  or  loss  of  wetland  func- 
tions such  as  flood  control,  nutrient  retention,  and 
wildlife  habitat. 

This  report  focuses  on  the  Gallatin  Valley  and  sur- 
rounding area,  typical  of  many  rapidly  growing  re- 
gions in  the  West  with  increasing  land  conversion, 
subdivision,  and  residential  development.  Our 
objectives  were  to  quantify  changes  in  wetland  ex- 
tent and  fiinction  in  our  study  area  and  to  estimate 
cumulative  change  in  wetland  functions. 

The  project  required  us  to  produce  new  digital 
wetland  maps  at  a  1:12,000  scale,  using  2005  aerial 
imagery  at  1 -meter  resolution.  This  was  done  as 
part  of  National  Wetland  Inventory  (NWI)  updat- 
ing, following  current  federal  standards. 

To  analyze  wetland  change,  we  compared  ran- 
domly selected  wetlands  from  the  original  NWI, 
completed  in  1984  and  1988,  with  the  new  NWI 
mapping  created  for  this  project.  We  randomly 
selected  25%  of  the  one-square  mile  Public  Land 
Survey  System  sections  in  each  subwatershed  in 
the  study  area  using  a  spatially  balanced  random 
sampling  approach.  Within  the  sampled  area,  we 
compared  each  wetland  polygon  in  the  old  map- 
ping to  the  corresponding  wetland  polygon  in  the 
new  mapping,  and  we  assigned  a  source  of  change 
to  each  mapped  wetland.  To  assess  the  func- 
tions associated  with  each  wetland,  we  analyzed 
the  landscape  position,  landform,  waterbody,  and 
water  flow  paths  for  each  wetland.  We  assigned 
hydrogeomorphic  (HGM)  modifiers  to  all  wetland 
polygons  in  both  the  old  and  new  wetland  mapping. 
These  HGM  attributes  were  combined  with  the 
NWI  classification  attributes  to  yield  a  combination 
that  could  then  be  ranked  on  a  performance  scale  of 
1  (low),  2  (moderate),  and  3  (high)  for  each  often 


wetland  functions  (water  storage,  streamflow  main- 
tenance, groundwater  recharge,  nutrient  cycling, 
sediment  retention,  shoreline  stabilization,  native 
plant  community  maintenance,  tertestrial  habitat, 
aquatic  habitat,  and  conservation  of  wetland  bio- 
diversity). We  used  this  performance  ranking  as 
a  weighting  factor  and  muhiplied  this  weighting 
factor  by  wetland  area  to  calculate  functional  units 
for  each  wetland  function.  We  also  completed  a 
wetland  landscape  profile  for  each  sixth  code  hy- 
drologic  unit  that  provides  a  broad  landscape  char- 
acterization of  wetlands  across  the  project  area. 

We  digitized  56,822  acres  (22,995  hectares)  of  wet- 
lands and  28,210  acres  (1 1,416  hectares)  of  riparian 
habitat  within  the  change  detection  analysis  area. 
Palustrine  emergent  wetlands  covered  the  greatest 
area  with  over  28,380  acres  ( 1 1 ,485  hectares).  The 
majority  of  wetland  and  riparian  habitats  (57,358 
acres;  23,212  hectares)  occur  on  private  lands 
within  the  analysis  area.  Overall,  we  observed  an 
increase  of  4,22 1  mapped  wetland  acres  ( 1 ,708 
hectares)  between  1980's  and  2005  within  the  study 
area. 

Wetlands  associated  with  lotic  features  comprised 
the  largest  hydrogeomorphic  type  in  the  Gallatin 
project  area,  totaling  39,454  acres  (15,967  hect- 
ares). Wetlands  associated  with  deepwater  and 
associated  lentic  features  covered  824  acres  (333 
hectares),  and  terrene  wetlands  totaled  3,256 
acres  ( 1 ,3 1 8  hectares).  Comparison  of  wetland 
functional  performance  capacities  throughout  the 
Gallatin  project  area  showed  an  overall  73.5%  gain 
in  hydrologic  functions  that  include  water  stor- 
age, streamflow  maintenance,  and  groundwater 
recharge.  However,  we  mapped  over  five  times 
more  acres  of  lotic  wetlands  using  higher  resolu- 
tion 2005  imagery,  which  contributed  to  this  appar- 
ent gain  in  hydrologic  fiinction.  Biogeochemical 
functions  incorporating  nutrient  cycling,  sediment 
retention,  and  shoreline  stabilization  showed  an 
overall  increase  of  24%.  Functions  associated  with 
terrestrial  and  native  plant  communities  showed  a 
combined  decline  of  7.3%,  whereas  aquatic  habitat 
and  conservation  of  wetland  biodiversity  showed  a 
combined  increase  of  9.7%. 


Ill 


Our  analysis  shows  an  overall  increase  in  wetland 
area  between  mapping  completed  in  the  1980"s  and 
new  mapping  from  2005  aerial  imagery.  Many 
wetlands  mapped  as  palustrine  scrub  shrub  in  the 
1980"s  are  now  palustrine  emergent  wetlands.  Ex- 
amination of  the  aerial  imagery  revealed  that  much 
of  this  change  is  attributable  to  agricultural  changes 
(e.g.,  livestock  grazing,  stream  dewatering,  and 
conversion  to  hay  pasture).  Additionally,  some 
scrub  shrub  historically  mapped  as  wetland  was 
mapped  as  riparian  scrub  shrub  in  the  new  map- 
ping. 

Although  our  estimates  indicate  that  few  actual 
wetland  acres  were  lost  between  the  1980's  and 
2005  mapping  efforts  across  the  entire  project  area, 
concentrated  wetland  losses  occurred  in  a  few 
areas.  In  particular,  much  of  the  wetland  change 
in  the  areas  immediately  around  Bozeman  was  at- 
tributable to  urban  and  rural  development.  This 
area  has  seen  rapid  growth  and  much  of  the  valley 
bottom  locations  along  the  East  Gallatin  River  have 
been  subdivided.  The  wetland  landscape  profil- 
ing also  revealed  that  the  areas  around  Bozeman 
contain  wetlands  with  the  potential  for  high  per- 
fonnance  of  several  wetland  functions,  including 
groundwater  recharge,  streamflow  maintenance, 
water  storage,  sediment  retention,  and  terrestrial 
habitat.  Continued  impacts  to  wetlands  in  these 
important  areas  will  reduce  the  ability  of  wetlands 
to  perfomi  these  functions,  potentially  resulting  in 
ecological  and  economic  losses  in  these  areas. 


It  is  also  important  to  note  that  differences  between 
the  scales  of  the  imagery  used  in  the  historic  and 
new  mapping  products  make  accurate  quantifica- 
tion of  wetland  change  problematic.  The  historic 
wetland  mapping  was  digitized  at  a  1 :58.000 
scale  and  exhibits  considerably  more  spatial  error 
than  the  current  mapping  that  was  digitized  at  a 
1:12.000  scale. 

Factors  such  as  photo  quality,  scale,  and  environ- 
mental conditions  at  the  time  of  photo  acquisition 
can  also  affect  mapping  accuracy.  Digital  wetland 
maps  are  static  and  may  not  reflect  the  dynamic 
nature  of  wetlands  subject  to  drastic  annual  and 
seasonal  fluctuations  in  size  and  distribution.  We 
also  emphasize  that  the  functional  capacity  ratings 
assigned  to  wetlands  in  this  project  are  only  poten- 
tial capacities.  Data  on  actual  functional  capacity 
would  require  extensive  field  work  and  assessment. 

This  analysis  should  be  considered  a  preliminary 
assessment  of  changes  in  the  Gallatin  Valley  wet- 
lands and  wetland  functional  capacity.  Data  from 
this  analysis  can  provide  ver\'  effective  conserva- 
tion tools  to  identify  areas  with  the  potential  to  per- 
form wetland  functions  most  effectively,  allowing 
natural  resource  managers  and  other  stakeholders 
to  focus  or  prioritize  conservation  and  restoration 
efforts. 
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Introduction 


Digital  wetland  mapping  provides  important  infor- 
mation on  the  t\'pe,  location,  and  extent  of  wetlands 
within  a  given  region.  Comparing  historic  wetland 
mapping  with  updated  mapping  provides  an  oppor- 
tunit\  to  examine  potential  changes  in  wetland  den- 
sity and  distribution  due  to  both  natural  and  anthro- 
pogenic causes.  The  National  Wetland  Inventory 
(NWl)  of  the  U.S.  Fish  and  Wildlife  Service  has 
completed  several  reports  on  the  status  and  trends 
of  the  nation's  wetlands  (Frayer  et  al.  1983,  Dahl 
and  Johnson  1991,  Dahl  2000,  2006),  underscoring 
the  cumulative  change  in  wetland  area  on  a  na- 
tional scale.  However,  documentation  of  wetland 
change  at  a  local  or  regional  scale  is  also  necessary, 
particularly  for  areas  undergoing  rapid  land  use 
changes.  In  addition  to  documenting  changes  in 
wetland  area,  addressing  the  corresponding  change 
or  loss  of  wetland  functions  such  as  flood  control, 
nutrient  retention,  and  wildlife  habitat  is  also  im- 
portant, as  wetlands  are  valued  not  for  the  area  they 
cover  but  for  the  ecological  functions  they  provide. 

Certain  wetland  functions  become  increasingly 
important  in  areas  influenced  by  human  develop- 
ment (Mitsch  and  Gosselink  2000a).  Yet,  wetland 
functional  perfomiance  can  vary  both  annually  and 
seasonally.  Moreover,  not  all  wetlands  perform 
all  functions.  For  example,  a  wetland  located  near 
a  headwater  stream  performs  different  functions 
than  a  wetland  at  the  stream's  mouth  (Mitsch  and 
Gosselink  2000b).  Thus,  each  wetland  is  unique  in 
both  its  ecology  and  functionality. 


wetlands.  The  hydrogeomorphic  classification 
(HGM)  was  developed  to  classify  wetlands  ac- 
cording to  their  geomorphic  position,  water  source, 
and  hydrodynamics  (Brinson  1993).  These  three 
factors  affect  how  a  wetland  functions:  and  wet- 
lands with  similar  characteristics  typically  function 
similarly.  The  addition  of  HGM  descriptors  to 
each  wetland  polygon  in  an  NWl  digital  database 
enhances  the  information  about  each  wetland.  Wet- 
land type  can  be  linked  to  wetland  function,  pro- 
viding a  simple  way  to  identify  wetlands  with  the 
potential  to  perform  various  ecological  functions 
(Tiner2005). 

In  an  assessment  context.  HGM  applications  are 
largely  site  specific,  based  on  regional  reference 
conditions,  and  require  substantial  field  investiga- 
tions. However,  recent  studies  in  the  northeastern 
U.S.  (Tiner  2005)  and  Colorado  (Johnson  2005) 
have  illustrated  the  application  of  HGM  as  a  tool 
for  initial,  coarse-scale  characterization  and  assess- 
ment of  wetland  functions  through  GIS-based  land- 
scape analyses.  This  tool  can  be  used  to  establish 
a  baseline  for  wetland  functionality,  and  to  create 
a  landscape-level  wetland  functional  profile.  This 
baseline  can  be  compared  with  current  conditions 
to  provide  information  on  the  status  and  trends  of 
wetland  functions,  allowing  for  the  identification  of 
potentially  ecologically  significant  wetlands,  pri- 
oritization of  restoration  and  protection  efforts,  and 
detection  of  wetlands  most  effective  at  performing 
particular  wetland  functions. 


Although  individual  wetlands  perform  unique  func- 
tions, they  do  not  act  in  isolation  but  are  instead 
influenced  by  the  surrounding  landscape.  Wetlands 
also  influence  the  function  and  condition  of  entire 
watersheds.  Mapping  wetland  resources  provides 
information  on  the  density  and  distribution  of 
wetlands  across  these  large  landscapes.  NWl  data 
are  the  most  widely  available,  comprehensive,  and 
standardized  wetland  mapping  products  available. 
The  NWl  classification  system  (Cowardin  et  al. 
1979)  also  provides  a  means  of  describing  wetland 
characteristics  such  as  water  regime,  vegetation, 
and  substrate.  However.  NWl  descriptors  provide 
little  information  on  the  functionality  of  individual 


Estimating  change  in  wetland  extent  and  distribu- 
tion as  well  as  change  in  wetland  function  over 
time  is  particularly  valuable  for  areas  experiencing 
rapid  land  use  changes.  In  Montana  and  much  of 
the  western  U.S..  wetlands  are  typically  concen- 
trated in  broad  river  valleys  and  riparian  areas.  The 
attractive  amenities  offered  by  such  areas  often 
lead  to  increased  development  and  land  conversion. 
These  pressures  can  have  negative  impacts  on  wet- 
land functions  including  reduced  groundwater  re- 
charge, increased  sediment  and  nutrient  inputs,  in- 
creased invasion  of  nonnative  species,  and  reduced 
wildlife  habitat. 


Increasing  concern  over  the  potential  loss  of  wet- 
lands and  their  associated  functions  due  to  rapid 
development  in  large  river  valleys  has  prompted 
examination  of  how  wetland  extent,  distribution, 
and  function  have  changed  in  these  areas.  This 
report  is  the  final  of  three  reports  focusing  on  wet- 
land change  and  ecological  function  in  rapidly  de- 
veloping valleys  of  Montana  (Kudray  and  Schemm 
2008,  Newlon  and  Bums  2009).  This  report  fo- 
cuses on  the  Gallatin  Valley  and  surrounding  area. 


which  is  typical  of  many  rapidly  growing  regions 
in  the  West  experiencing  increased  land  conver- 
sion, subdivision,  and  residential  development. 
Our  objectives  were  to  quantify  changes  in  wet- 
land extent  and  function  in  our  study  area  and  to 
estimate  cumulative  change  in  wetland  functions. 
In  conjunction  with  new  NWI  maps  produced  by 
the  MTNHP.  this  report  provides  a  comprehensive 
characterization  of  wetlands  in  the  study  area. 


Study  Area 


The  Gallatin  Valley  study  area  includes  portions 
of  Broadwater,  Jefferson,  Madison,  and  Gallatin 
counties  in  the  Middle  Rockies  Level  III  Ecoregion 
(Omemik  1987).  Elevation  ranges  from  4,000  feet 
( 1,2 19  meters)  in  the  valley  bottom  near  Logan 
to  9,012  feet  (2,747  meters)  in  the  Tobacco  Root 
Mountains.  The  area  mapped  is  over  967,000  acres 
(391.000  hectares)  and  includes  portions  of  70  sub- 
watersheds  (sixth  code  hydrologic  units  or  HUCs) 
covering  29  U.S.  Geological  Survey  1 :24,000  topo- 
graphic quadrangles  (Figure  1 ).  Outside  of  urban 
areas,  land  use  is  predominantly  wildlife  habitat 
and  recreation.  Livestock  grazing  and  farming  oc- 
cur in  valley  locations.  Approximately  86%  of  the 
analysis  area  is  privately  owned  (Figure  2).  The 
largest  city  is  Bozeman  (Gallatin  County). 

Climate 

Regional  climate  is  influenced  by  continental 
air  masses  with  cold,  sunny,  and  dry  winters  and 
warm,  dry  summers.  In  the  valley  bottoms  around 
Bozeman,  average  yearly  maximum  temperature  is 
55.2°F  (12.9°C)  and  average  yearly  minimuin  tem- 
perature is  3 1 .2°F  (-0.4°C).  The  warmest  months 
are  July  and  August.  Average  annual  precipitation 
is  18.4  inches  (467  mm).  A  third  of  the  yearly 
precipitation  occurs  during  April,  May,  and  June. 
Average  total  snowfall  is  85  inches  (2,162  mm) 
(Western  Regional  Climate  Center  2009). 

Geology 

The  area  is  characterized  by  wide  intermontane  val- 
leys and  several  mountain  ranges.  Geology  of  the 
project  area  is  complex.  The  majority  of  the  proj- 
ect area  is  comprised  of  large  river  valleys  primar- 
ily fornied  in  Tertiary  sediments  and  Quaternary 
alluvial  deposits  (Nesser  et  al.  1997).  Mountain 
ranges  in  the  area  were  largely  fornied  in  mixed 
sedimentary  or  metasedimentary  bedrock,  with  the 


exception  of  the  Gallatin  Mountains  in  the  south- 
eastern portion  of  the  project  area,  which  were 
formed  in  Tertiary  volcanic  flows  and  associated 
pyroclastic  deposits  consisting  of  rhyolite,  basalt, 
and  andesite. 

Hydrology 

The  study  area  encompasses  portions  of  the  follow- 
ing watersheds:  the  Jefferson  River,  the  Madison 
River,  the  Gallatin  River,  and  the  East  Gallatin 
River  (Figure  3).  Primary  water  uses  are  irrigation, 
livestock,  and  domestic  wells. 

Vegetation 

Natural  vegetation  communities  are  the  most  com- 
mon land  cover  type  in  the  project  area  (Table  1 ). 
Although,  pasture  and  cultivated  cropland  comprise 
nearly  25%  of  the  area.  Land  cover  and  land  use 
data  for  each  sixth  code  hydrologic  unit  are  shown 
in  Appendix  G.  Montane  portions  of  the  project 
area  are  dominated  by  coniferous  forests.  Higher 
elevation  forests  are  comprised  of  whitebark  pine 
{Pinus  alhicaulis),  subalpine  fir  {Abies  lasiocarpa), 
Engelmann  spruce  (Picea  engelmannii),  and  lodge- 
pole  pine  {Pinus  contorta).  Middle  and  lower 
elevation  forests  are  dominated  by  Douglas-fir 
{Pseudotsuga  meiiziesii)  and  lodgepole  pine.  Lim- 
ber pine  {Pinus  flexilis)  and  Rocky  Mountain  ju- 
niper {Juniperus  scopulorum)  are  locally  common 
on  calcareous  substrates.  Quaking  aspen  {Popuhis 
tremuloides)  occurs  in  stands  throughout  the  proj- 
ect area  in  sites  where  adequate  soil  moisture  is 
available. 

Foothills  and  valleys  are  largely  grassland  and 
shrubland.  Mountain  big  sagebrush  {Artemisia 
tridentata  ssp.  vaseyana)  is  the  dominant  shrub, 
although  stands  of  threetip  sagebrush  {A.  tripar- 
tita) and  basin  big  sagebrush  {A.  tridentata  ssp. 


Table  1.  Land  cover  and  land  use  area  (in  acres)  for  the  Gallatin  project  area  from  the  2009  Mon- 
tana Land  Cover  theme. 
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tridentata)  are  locally  common.  Alkaline  or  saline 
sites  may  support  stands  of  greasewood  (Sarco- 
batus  vermiciilatiis).  Bunchgrasses  dominate  the 
herbaceous  layer,  including  bluebunch  wheatgrass 
{Pseiidoroegneria  spicata),  Idaho  fescue  (Festuca 
idahoensis),  needle-and-thread  (Hesperostipa  co- 
niata),  Sandberg's  bluegrass  {Poa  secinida),  and 
blue  grama  {Boutehua  gracilis). 

A  previous  study  identified  high  quality  wetlands  in 
this  area  (Jones  2004).  The  following  descriptions 
of  wetland  and  riparian  types  in  the  Gallatin  project 
area  are  based  on  observations  from  that  report. 

Riparian  and  Palustrine  Forested 
Wetlands 

The  large  river  valleys  consist  of  riparian  and 
palustrine  forests  dominated  by  black  cottonwood 
{Popu/us  balsanufera)  (Figure  4).  Cottonwood 
stands  vary  from  open  woodlands  to  closed  canopy 
forests.  Some  areas  along  the  lower  reaches  of  the 
Jefferson.  Madison,  and  Gallatin  Rivers  still  have 
cottonwood  recruitment,  although  the  shrub  layer 
has  been  lost  in  many  areas.  On  those  floodplains 
with  a  well-developed  shrub  layer,  yellow  willow 
(Salix  lutea),  redosier  dogwood  (Conms  sericea), 
common  chokecherry  {Primus  virginiana).  west- 
em  snowberry  {Symphoricarpos  occidentalis).  and 
Wood's  rose  (Rosa  woodsii)  are  found.  Degraded 
stands  have  little  to  no  shrub  development  or  a 
shnib  layer  dominated  by  western  snowberry  and 
Wood's  rose.  The  herbaceous  layer  is  largely  dom- 
inated by  non-native  pasture  grasses,  particularly 


Kentucky  bluegrass  {Poa  pratensis),  smooth  brome 
(Brormis  inermis),  and  common  timothy  (Phleum 
pratense).  Higher  elevation  streams  are  dominated 
by  Engelmann  spruce,  subalpine  fir,  and  quaking 
aspen. 

Riparian  and  Palustrine  Scrub  Shrub 
Wetlands 

Riparian  and  wetland  scrub  shrub  vegetation  oc- 
curs on  terraces  along  the  floodplains  of  both  low 
and  high  gradient  streams  and  rivers  (Figure  5).  It 
can  also  be  found  around  beaver  ponds  and  on  the 
fringes  of  fens  and  lakes.  Common  willow  species 
at  higher  elevations  include  planeleaf  willow  {Salix 
planifolia)  and  Wolf's  willow  (5'.  wolfii).  At  middle 
elevations,  Drummond's  willow  {S.  dnimmondi- 
ana),  Geyer  willow  {S.  geyeriana),  and  Booth's 
willow  (5*.  boothii)  are  common.  At  lower  eleva- 
tions, yellow  willow  and  sandbar  willow  {S.  ex- 
igua)  occur  along  floodplains  and  on  active  gravel 
bars.  Other  common  shrub  species  include  moun- 
tain a\dQv{Alnus  incana),  chokecherry.  Wood's 
rose,  and  western  snowberry. 


Figure  4.  Cottonwood  forest  along  tlie  Gallatin  River 


Figure  5.  Palustrine  scrub  shrub  wetland  ui  the  Oatlalin  proj- 
ect area. 

Palustrine  Emergent  Wetlands 

Herbaceous  emergent  vegetation  in  the  project  area 
is  found  in  a  variety  of  settings  with  several  sedge 
species  dominant,  including  Carex  lasiocarpa,  C. 
utriculata,  C.  vesicaria,  and  C.  limosa  occurring 
on  the  wettest  sites.  Native  grasses  such  as  tufted 
hairgrass  {Deschampsia  cespitosa)  and  bluejoint 
reedgrass  {Calamagrostis  canadensis)  are  com- 


mon  on  wet  meadow  sites  (Figure  6).  Wet  meadow 
communities  iiave  been  greatly  impacted  by  live- 
stock and  other  agricultural  practices  are  largely 
dominated  by  non-native  grasses,  particularly  Ken- 
tucky bluegrass.  smooth  brome,  common  timothy, 
redtop  {Agrostis  gigantea),  and  reed  canarygrass 
(Phalaris  arundinacea). 


Figure  6.  Paliistrine  emergent  M-etlands  in  the  Gallatin  prajcci 
area. 

Aquatic  Bed 

Palustrine,  Lacustrine,  and  Riverine  aquatic  bed 
wetlands  occur  in  littoral  (<  2  meters  in  depth)  and 
limnetic  (>  2  meters  in  depth)  zones  of  ponds  and 
lakes  (Figure  7)  or  on  the  beds  of  slow-moving  pe- 
rennial streams.  Common  aquatic  plant  species  in- 
clude the  floating-leaved  yellow  pond  lily  {Nuphur 
polysepaluDi),  submerged  or  partially  submerged 
species  such  as  water  milfoil  (Myriophylliini  spp.). 
common  mare's-tail  (Hippiiris  vulgaris),  and  pond- 
weeds  (Potamogeton  spp.). 


Figire  7.  Palusii  nic  lujiuUic  hcd  n  liUuhI  hi  iIic  ( ml  latin  proj- 
ect area. 


Land  Use  History 

The  following  historical  infonnation  is  taken  from 
Smith  (1996).  The  Blackfeet  tribe  considered  the 
Gallatin  Valley  to  be  prime  hunting  grounds.  Later, 
other  tribes  including  the  Crow,  Flathead,  and  Sho- 
shone began  to  use  the  area,  as  did  early  European 
explorers,  fur  trappers,  and  settlers.  The  expedition 
led  by  Lewis  and  Clark  through  the  area  in  1805 
signaled  the  beginning  of  settlement  in  the  region. 
However,  the  discovery  of  gold  southwest  of  the 
Gallatin  Valley  in  the  mid-1860's  resulted  in  the 
establishment  of  the  first  permanent  settlements  in 
the  area. 

One  of  the  first  settlements  occurred  in  the  area  at 
present  day  Three  Forks,  where  agriculture  became 
an  important  land  use.  Crops  found  markets  in  the 
nearby  mining  communities  of  Virginia  City  and 
Helena.  Livestock  operations  had  an  important 
economic  impact  on  the  early  history  of  the  area 
as  well.  Settlements  formed  in  the  Gallatin  Valley 
around  1 864  at  the  site  of  present  day  Bozeman. 
The  Gallatin  Valley  became  a  transportation  hub 
with  the  arrival  of  the  Northern  Pacific  Railroad  in 
1 883.  Livestock  grazing  and  agriculture  remain 
important  land  uses  in  the  area. 

The  Gallatin  area  has  seen  rapid  population  growth 
over  the  past  century,  particularly  over  the  last  30 
years.  The  populations  of  all  four  counties  have 
increased  since  1980  (Figure  8).  Growth  has  been 
especially  rapid  in  Gallatin  County,  where  the 


population  more  than  doubled  since  1980  (U.S. 
Census  Bureau  2009).  Population  growth  in  ail 
counties  except  Madison  has  outpaced  that  of  both 
Montana  and  the  nation  (Headwaters  Economics 
2009a,  2009b.  2009c,  2009d).  Similarly,  density 
of  single-family  housing  has  increased  since  1980, 


particularly  in  Gallatin  County  (Figure  9).  A  com- 
parison of  the  number  of  water  rights  between  1980 
and  2008  illustrates  the  increased  demand  on  water 
resources  within  the  Gallatin  Valley,  particularly 
around  the  cities  of  Manhattan,  Belgrade,  and  Boz- 
eman  (Figure  1 0). 
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Figures.   Population  growth  from  1980-2008  in  Gallatin.  Madison.  Broadwater  and 
Jefferson  Counties.  Source:  U.S.  Census  Bureau  2009. 
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Figure  9.  Growth  in  single-family  housing  density  in  Gallatin.  Madison.  Broadwater, 
and  Jefferson  Counties  from  1980-2006.  Source:  Headwaters  Economics  2009. 
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Methods 


Wetland  and  Riparian  Mapping 

To  create  the  new  wetland  and  riparian  mapping, 
we  digitized  wetland  and  riparian  areas  using  Arc- 
GIS  9.3  (ESRI  2008).  We  used  1:12.000  scale  col- 
or-infrared aerial  photography  taken  in  2005  along 
with  several  ancillary  data  layers  (Table  2)  to  en- 
hance mapping  accuracy.  The  NWI  classification 
system  (Cowardin  et  al.  1979)  was  used  to  classify 
wetlands.  Mapping  procedures  followed  the  US- 
FWS  National  Standards  and  Quality  Components 
(USFWS  2004a)  and  the  Technical  Procedures  for 
Mapping  Wetland,  Deepwater,  and  Related  Habi- 
tats (USFWS  2004b).  Historic  wetland  mapping 
was  completed  for  the  project  area  by  the  National 
Wetlands  Inventory  using  1:58,000  scale  color-in- 
frared aerial  photography  in  1984  and  1988  (Figure 
11). 

The  NWI  classification  system  is  a  hierarchical 
approach  that  classifies  wetlands  into  systems, 
subsystems,  and  classes  (Appendix  A).  Three  sys- 
tems occur  in  Montana:  Palustrine  wetlands  (e.g., 
marshes,  fens,  potholes,  wet  meadows,  willow 
carrs.  forested  swamps);  Lacustrine  wetlands  (e.g., 
associated  with  waterbodies  at  least  6.6  feet  [2  me- 
ters] in  depth  and/or  at  least  20  acres  [8  hectares]  in 
size);  and  Riverine  wetlands  (e.g.,  occurring  within 
the  active  channel  of  a  stream  or  river).  Lacustrine 
and  Riverine  svstems  are  further  classified  into 


subsystems.  Lacustrine  subsystems  include  lim- 
netic (>  2  meters  in  depth)  and  littoral  (<  2  meters 
in  depth)  zones.  Riverine  subsystems  that  occur  in 
Montana  include  Lower  Perennial,  Upper  Peren- 
nial, and  Intemiittent.  Classes  describe  the  sub- 
strate for  Lacustrine  and  Riverine  systems,  and  the 
vegetation  life  forni  (e.g.,  emergent,  scrub-shrub, 
forested)  for  Palustrine  wetlands.  Hydrologic  mod- 
ifiers are  applied  to  each  wetland  to  describe  the 
water  regime  along  a  continuum  from  wettest  (per- 
manently flooded)  to  driest  (temporarily  flooded). 
Additional  modifiers  are  added  to  wetlands  altered 
by  anthropogenic  (e.g.,  excavated,  impounded)  and 
natural  (e.g.,  beaver)  influences. 

Riparian  habitats  have  not  been  previously  mapped 
in  Montana.  We  used  the  U.S.  Fish  and  Wildlife 
Service  Western  Riparian  Classification  System 
(USFWS  1997;  Appendix  B)  to  classify  riparian 
areas.  Riparian  areas  typically  lack  the  amount  or 
duration  of  water  present  in  wetlands,  yet  are  more 
mesic  and  vegetation  is  more  vigorous  than  adja- 
cent uplands.  Woody  vegetation  is  the  predomi- 
nant fype  of  riparian  habitat  mapped  (e.g..  riparian 
shrubland  and  riparian  forest),  although  herbaceous 
riparian  emergent  vegetation  may  be  mapped  if 
the  imagery  allows  for  identification  of  this  fype. 
The  riparian  classification  system  is  similar  to  the 
Cowardin  classification  system  in  its  hierarchical 
approach.  The  riparian  classification  system  has 


Table  2.  Ancillan'  spatial  data  layers  used  to  digitize  and  classify  wetlands  and  riparian  areas. 


Data  Layer 

Data  Source 

Year 

1 -meter  resolution  color- infrared  digital 
orthophotography 

National  Agricultural  Imagery  Program 

2005 

1 -meter  resolution  true  color 

National  Agricultural  Imagery  Program 

2005 

1 -meter  black  and  white  digital  orthophoto  quarter 
quadrangles 

U.S.  Geological  Survey 

1990 

10-meter  digital  elevation  model 

U.S.  Geological  Survey 

Various 

1 :24,000  topographic  quadrangle  digital  raster 
graphic 

U.S.  Geological  Survey 

Various 

1:24,000  high  resolution  National  Hydrography 
Dataset 

U.S.  Geological  Survey 

Various 

Original  National  Wetlands  Inventory  Wetland 
Mapping 

National  Wetlands  Inventory,  U.S.  Fish 
and  Wildlife  Service 

1981 
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a  single  system  (Riparian),  two  subsystems  (Lotic 
and  Lentic).  and  three  classes  (Forested,  Scrub- 
shrub,  and  Emergent). 

Wetland  Change  Detection 
Analysis 

Our  sampling  units  were  one-square  mile  Public 
Land  Survey  System  (PLSS)  sections.  We  ran- 
domly selected  25%  of  the  sections  in  each  sixth 
code  hydrologic  unit  in  the  study  area  using  a  spa- 
tially balanced  random  sampling  approach.  We 
used  a  Generalized  Random  Tessellation  Stratified 
(GRTS)  sampling  design  to  select  PLSS  sections 
(Stevens  1997,  Stevens  and  Olsen  1999,  Stevens 
and  Olsen  2004),  using  the  "spsurvey"  package 
for  Spatial  Survey  Design  and  Analysis  (Kincaid 
et  al.  2008)  in  R  2.8.1  (R  Core  Development  Team 
2008).  Historic  wetland  mapping  that  occurred 
within  these  selected  PLSS  sections  was  used  in  the 
change  detection  analysis. 

Improved  image  quality  and  resolution  associated 
with  the  2005  imagery  has  enhanced  our  ability 
to  map  wetlands  more  accurately  by  delineating 
wetland  boundaries  more  precisely  or  detecting 
wetlands  that  may  have  been  excluded  from  previ- 
ous mapping  eflforts.  We  made  every  effort  to  cor- 
rectly recognize  and  attribute  the  sources  of  wet- 
land change  by  visually  inspecting  each  wetland 


polygon  in  both  the  old  and  new  wetland  mapping 
within  the  sampled  sections.  Within  the  sampled 
area,  we  compared  each  wetland  polygon  in  the 
old  mapping  to  the  corresponding  wetland  polygon 
in  the  new  mapping,  and  we  assigned  a  source  of 
change  to  each  polygon  based  on  the  USFWS  Sta- 
tus and  Trends  Procedures  (2005)  with  some  modi- 
fications (Table  3).  For  example,  if  a  wetland  was 
mapped  as  palustrine  emergent  in  the  old  mapping 
but  had  been  excavated  to  create  a  pond,  the  poly- 
gon was  attributed  as  a  "New  Pond."  If  the  new 
wetland  polygon  was  classified  the  same  as  the  old 
wetland  polygon,  then  the  polygon  was  attributed 
with  "No  Change."  We  also  included  a  "Mapped 
as  Riparian"  category  because  riparian  areas  were 
not  previously  mapped.  Wetland  changes  attrib- 
utable to  interpreter  differences  or  differences  in 
image  quality  were  attributed  with  "Interpreter  Dif- 
ference" and  removed  from  the  analysis.  Addition- 
ally, we  visually  inspected  all  wetlands  associated 
with  beaver  ponds  and  all  created  wetlands  (e.g., 
excavated  ponds). 

Wetland  water  regimes  can  change  annually,  so 
we  analyzed  wetland  change  to  the  class  level 
only.  We  estimated  the  difference  in  wetland  acre- 
age between  the  old  and  new  mapping  using  the 
"spsurvey"  package  for  Spatial  Survey  Design  and 
Analysis  (Kincaid  et  al.  2008)  in  R  2.8. 1  (R  Core 
Development  Team  2008). 


Table  i.   Change  categories  assigned  to  sampled  wetland  polygons. 


Code 

Type 

Description 

AG 

Agriculture 

Land  used  for  the  production  of  food/fiber.  Often  evidenced  by 
geometric  patterns. 

RD 

Rural  Development 

Non-intensive  land  use  and  sparse  building  density. 

UD 

Urban  Development 

Land  is  predominantly  covered  by  structures  (high  building  density). 

NC 

No  Change 

No  change  in  wetland  attribute. 

NA 

Natural 

Change  in  wetland  class  due  to  natural  changes  in  serai  stage  (e.g., 
PEM  succeeding  to  PSS). 

RE 

Resource  Extraction 

Oil  and  gas,  mining,  wind  farm. 

OT 

Other 

Areas  in  transition  (filled,  drained),  but  final  land  use  is  undetermined 
at  date  of  imagery. 

IN 

Interpreter  differences 

Wetland  change  due  to  interpreter  error,  differences  in  imagery,  etc. 

PD 

New  pond 

Pond  created  by  either  excavation  or  impoundment. 

RP 

Mapped  as  Riparian 

Wetland  mapped  as  riparian  in  new  mapping 
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Wetland  Functional  Performance 

To  assess  the  functions  associated  with  each  wet- 
land, we  analyzed  its  landscape  position,  land- 
form,  waterbody.  and  water  flow  path.  We  used 
an  HGM  attribute  coding  approach  developed  by 
Tiner  (2003.  2005)  and  modified  by  the  Montana 
Natural  Heritage  Program.  This  approach  has  been 
demonstrated  in  previous  wetland  analyses  (Vance 
et  al.  2006,  Kudray  and  Schemm  2008,  Vance 
2009)  but  has  been  further  refined  for  this  study 
(Appendices  C  and  D).  We  assigned  an  HGM  at- 
tribute code  to  all  wetland  polygons  in  both  the  old 
and  new  wetland  mapping  (Appendix  E).  These 
HGM  attributes  were  combined  with  the  NWI  clas- 
sification attributes  to  yield  a  combination  that  can 
then  be  ranked  on  a  performance  scale  of  1  (low). 
2  (moderate),  and  3  (high)  for  each  often  wetland 
functions  (water  storage,  streamflow  maintenance, 
groundwater  recharge,  nutrient  cycling,  sediment 
retention,  shoreline  stabilization,  native  plant  com- 
munity maintenance,  terrestrial  habitat,  aquatic 
habitat,  and  conservation  of  wetland  biodiversity) 
(Appendix  F).  We  used  this  performance  ranking 
as  a  weighting  factor  and  multiplied  this  weight- 
ing factor  by  wetland  area  to  calculate  functional 
units  for  each  wetland  fimction  (see  Tiner  2005). 


This  method  does  not  account  for  actual  wetland 
functional  perfonnance  but  rather  ranks  potential 
wetland  functional  performance  based  on  presumed 
functioning  condition.  Rationale  for  wetland  flinc- 
tional  performance  rankings  can  be  found  in  Vance 
et  al.  (2006)  and  Kudray  and  Schemm  (2008). 

Wetland  Landscape  Profiling 

A  wetland  landscape  profile  provides  a  broad  land- 
scape characterization  of  the  wetlands  within  a 
particular  area  at  the  subwatershed  level.  For  this 
project,  we  completed  analysis  at  the  sixth  code 
hydrologic  unit  level  using  the  new  mapping  based 
on  2005  imager},.  We  calculated  five  metrics  to 
produce  the  wetland  profile:  overall  wetland  acres; 
acres  of  isolated  wetlands  (defined  as  wetlands  not 
located  within  a  stream  or  on  its  floodplain  and 
not  connected  to  another  wetland);  acres  of  altered 
wetlands;  percent  of  wetlands  in  private  owner- 
ship; and  percent  of  wetlands  with  high  functional 
value  for  each  often  wetland  ecological  functions. 
For  each  function,  the  acres  of  wetland  assigned  a 
value  of  "3"  for  high  function  were  summed  and 
divided  by  the  total  number  of  acres  within  each 
subwatershed. 
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Results 


We  digitized  56,822  acres  (22,995  hectares)  of  wet- 
lands and  28.2 1 0  acres  ( 1 1 ,4 1 6  hectares)  of  riparian 
habitat  within  the  change  detection  analysis  area. 
Palustrine  emergent  wetlands  covered  the  greatest 
area  with  over  28,380  acres  ( 1 1 ,485  hectares)  (Fig- 
ure 12). 

The  majority  of  wetland  and  riparian  habitats 
(57,358  acres;  23,212  hectares)  occur  on  private 


lands  within  the  analysis  area.  Palustrine  emergent 
wetlands  occupied  the  largest  acreage  with  over 
25,749  acres  (10,420  hectares)  (Figure  13).  On 
public  lands,  over  6,126  acres  (2,479  hectares)  of 
wetlands  and  riparian  areas  were  mapped,  with 
lands  owned  by  The  Nature  Conservancy  having 
the  most  wetlands  (Figure  14). 


Figure  12.  National  Wetland  Inventory  (NWI)  wetlands  mapped  using  2005  color-in- 
frared aerial  pholograplty:  NIVI  wetland  types  with  acreage  too  low  to  present  graphi- 
cally are  included  in  the  graph  inset.  See  Appendix  A  for  NWI  codes. 


Lacusmnc  Ponds  Emergent      Scrub-shrub    Rivennc  and      Ripanan  Ripanan  Kiponao       Grand  Total 

WclJand  Wetland  Rj\cnnc         Emergent      Scrub-shrub         Forest 

Shore 
WeDand/Ripaiian  Type 

Figure  13.  Acreage  of  wetland  and  riparian  types  mapped  on  private  lands  within  the 
Gallatin  project  area. 


15 


Land  Onncnhip 

Figure  14.  Acreage  oj  all  wetland  and  riparian  types  combined  mapped  on  public 
lands  within  the  Gallatin  project  area 


Wetland  Change  Detection 
Analysis 

We  observed  an  increase  of  4.221  wetland  acres 
(1,708  hectares)  between  the  old  and  new  map- 
ping, although  estimates  were  highly  imprecise 
for  some  wetland  types  w  ith  very  few  polygons 
(Table  4).  Palustrine  emergent  wetlands  showed 
an  increase,  while  palustrine  scrub  shrub  wetlands 
showed  a  decline.  Many  of  the  wetlands  associated 
with  rivers  mapped  in  the  1980"s  were  mapped  as 
upper  perennial,  whereas  riverine  wetlands  in  the 
new  mapping  were  classified  as  lower  perennial. 
Additionally,  using  the  larger  scale  2005  imagery 
often  resulted  in  a  wetland  feature  increasing  size. 
Due  to  small  sample  sizes,  we  could  not  estimate 
wetland  change  separately  for  each  sixth  code  hy- 
drologic  unit. 

Within  the  wetlands  sampled,  698  acres  (282 
hectares)  mapped  as  wetland  in  the  historic  map- 
ping were  mapped  as  riparian  in  the  new  wetland 
mapping  (Figure  15).  Appro.ximately  10  acres  (4 
hectares)  of  beaver-created  wetlands  have  been  lost 
since  the  historic  mapping  was  completed  (Figure 
16).  Diked/impounded  palustrine  wetlands  largely 
showed  a  decline  across  the  project  area:  however, 
excavated  wetlands  showed  an  increase,  particu- 
larly in  the  form  of  canals  and  ditches. 


Wetland  Functional  Performance 

Wetlands  associated  with  lotic  features  comprised 
the  largest  hydrogeomorphic  type  in  the  Gallatin 
project  area,  totaling  39,454  acres  (15,967  hectares) 
(Table  5).  Wetlands  associated  with  deepwater  and 
associated  lentic  features  covered  824  acres  (333 
hectares),  and  terrene  wetlands  totaled  3,256  acres 
(1.318  hectares). 

Comparison  of  wetland  functional  perfonnance 
capacities  throughout  the  Gallatin  project  area 
showed  an  overall  73.5%  gain  in  hydrologic  func- 
tions that  include  water  storage,  streamflow  main- 
tenance, and  groundwater  recharge.  However,  we 
mapped  over  five  times  more  acres  of  lotic  wet- 
lands using  the  larger  scale  2005  imagery,  which 
largely  contributed  to  this  change  in  hydrologic 
function.  Biogeochemical  functions  incorporating 
nutrient  cycling,  sediment  retention,  and  shoreline 
stabilization  showed  an  overall  increase  of  24%. 
Functions  associated  with  terrestrial  and  native 
plant  community  maintenance  showed  a  combined 
decline  of  7.3%.  whereas  aquatic  habitat  and  con- 
servation of  wetland  biodiversity  showed  a  com- 
bined increase  of  9.7%  (Table  6). 

Wetland  Landscape  Profiling 

The  wetland  landscape  profiling  analysis  revealed 
that  the  confluence  of  the  Madison,  Jefferson,  and 
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Table  4.   Changes  in  wetland  area  within  the  Gallatin  project  area  (in  acres)  generalized  by  NWI  wetland  type  between  the  1980  's 
and  2005. 


198trs                                                    1 

NW 1  Class 

#of 

mapped 

polygons 

included 

in  analysis 

Estimated 
Acres 

Standard 
Error 

95%  Confidence 
Limits 

LI  LIB 

3 

2*^4  9 

79  7 

138  8 

451  0  1 

l:ab 

none  mapped                                          | 

l:us 

4 

180 

8  0 

2  3 

337 

PABPUB 

175 

631  4 

87  4 

460  2 

802  7 

PFM 

375 

^,502  9 

1,820  4 

5,934  9 

13.070  8 

PFO 

20 

117  5 

33  5 

51  8 

183  3 

PSS 

213 

3,160  9 

437  2 

2,313  0 

4,026  8 

PUS 

4 

10  7 

62 

-1  5 

23  0 

R2UB 

S-' 

2,273  1 

405  4 

1,478  5 

3,067  6 

R2US 

Q6 

1.129  8 

199  1 

739  6 

1,520  0 

R3UB 

3  b 

1,503  1 

225  6 

1.060  9 

1,945  4 

R.UIS 

53 

513  1 

96  8 

323  3 

7029 

R-4SB 

T 

1  0 

0(1 

0  9 

1  0 

All  Wetlands 

1033 

19,1655 

3,399  5 

12,502.6 

25,828.3 

2005                                                     1 

#of 

mapped 

polygons 

included 

in  analysis 

Estimated 
.Acres 

Standard 
Error 

95%  Confidence 
Limits 

Estimated 
Wetland 
Change 

5 

341  7 

132  2 

82  5 

600  9 

468 

3 

30  5 

13  3 

44 

566 

— 

2 

19 

02 

1  5 

24 

-16  1 

126 

538  1 

87  6 

366  3 

709  8 

-93  4 

280 

13,829  4 

2,733  9 

8,471  0 

19,187  7 

4,326  5 

none  mapped 

— 

122 

1,701  9 

290  9 

1,131  9 

2,272  0 

■1,468  0 

15 

12  0 

37 

47 

192 

1  2 

98 

5,102  4 

568  1 

3,989  0 

6.215  8 

2,829  3 

none  mapped 

— 

5 

786 

26  6 

265 

130.8 

-1,4245 

none  mapped 

— 

4 

20  2 

129 

-5  1 

45  6 

19  3 

660 

21,656.6 

3,869,5 

14,0726 

29,240.7 

4,221,1 

Agnculture  liilcr]>rrtet 

DilTercncc 


Pond  Rural  DCTelopmLtii    MuppeJ  as  Rjpannn    Urban  Developmenl 


T>'pe  of  Change 

Figure  15.   Types  of  wetland  changes  (in  acres)  between  historic  and  new  wetland 
mapping  in  the  Gallatin  project  area. 
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■  Diked^mpounded 
E3  E:xcavaled 
D  Beaver 


NWI  Code 

Figure  16.  Changes  in  naturally  modified  (beaver)  and  created  wetlands  (in  acres) 
between  historic  and  new  wetland  mapping  in  the  Gallatin  project  area. 


Table  5.  ffydrogeomorphic  types  and  area  (in  acres)  mapped  within  the  Gallatin  project  area  for  both  historic  (left)  and  tww 
wetland  mapping  (right I. 


Landscape  Position/ 
Waterbody  Type 

Landform 

Water 
Flow  Path 

Area 

(Acres) 

Deepwalcr 

1  lironglinow 

474,4 

River 

3635.9 

Stream 

858.1 

Lentic 

Basin 

Ihrougliflou 

37.4 

Fringe 

Bidirectional 

II  1.3 

Ihroughflovv 

2.3 

Island 

Bidirectional 

0.5 

Lotic  River 

Basin 

Ihroiigliflow 

1 73. 1 

Fioodplain 

Throughflovv 

3008.4 

Fringe 

Ihroughflovv 

1060.9 

Island 

Ihroughflow 

473.7 

Lolic  Stream 

Basin 

Throughflow 

6916.6 

Fringe 

Ihroughflow 

1614.9 

Island 

Ihroughflow 

81.8 

Terrene 

Basin 

Isolated 

472.3 

Complex 

479.1 

Slope 

Isolated 

285.0 

Complex 

I0I.2 

Landscape  Position/ 
Waterbody  Type 

Landform 

Water 
Flow  Path 

Area 

(Acres) 

Deepwater 

Ihroughflow 

529. 1 

River 

4.451.5 

Stream 

935.7 

Ixntic 

Basin 

Ihroughflow 

116.8 

Fringe 

Bidirectional 

169.9 

Ihroughflow 

7.7 

Lotic  River 

Basin 

Ihroughflow 

483.8 

Fioodplain 

Ihroughflow 

15.516.9 

Fringe 

Ihroughflow 

1,558.6 

Lotic  Stream 

Basin 

Ihroughflow 

16.120.5 

Fringe 

Ihroughflow 

387.5 

lerrenc 

Basin 

Isolated 

613.0 

Complex 

991.8 

Slope 

Isolated 

814.3 

Complex 

836.5 
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Table  6.  Predicted  change  in  wetland  functional  capacity  between  historic  and  new  wetland  mapping  for  the  Gallatin  project 
area. 


Function 

1980's 

Functional 

Units 

2005 

Functional 

Units 

Predicted  Percent 

Change  in  Functional 

Capacity 

Hydrology 

Water  Storage 

11.606.7 

12,998.9 

12.0 

Streamflow  Maintenance 

11,894.7 

14,577.3 

22.6 

Groundwater  Recharge 

14,144.1 

19,652.1 

38.9 

Biogeochemical 

Nutrient  Cycling 

12,079.7 

12,576.0 

4.1 

Sediment  Retention 

13,339.5 

15,204.2 

14.0 

Shoreline  Stabilization 

7,546.9 

7,995.4 

5.9 

Habitat 

Native  Plant  Community  Maintenance 

8,693.3 

8,269.7 

-4.9 

Terrestrial  Habitat 

9,311.3 

9,083.6 

-2.4 

Aquatic  Habitat 

11,837.3 

12,717.3 

7.4 

Conservation  of  Wetland  Biodiversity 

7,621.6 

7,796.0 

2.3 

Gallatin  Rivers  near  Three  Forks  has  the  highest 
density  of  wetland  acres  followed  by  the  area  along 
the  Gallatin  River  (Figure  1 7).  The  area  southeast 
of  Harrison  has  the  largest  acreage  of  altered  wet- 
lands, represented  largely  by  Willow  Creek  Reser- 
voir (Figure  18).  Lands  to  the  south  of  Belgrade 
and  west  of  Bozeman  have  the  highest  acreage 
of  isolated  wetlands  along  with  the  west  slope  of 
the  Bridger  Mountains  (Figure  i  9).  The  wetland 
landscape  profile  also  indicated  that  the  majority  of 
wetlands  mapped  within  the  project  area  are  locat- 
ed on  private  lands,  particularly  in  the  western  por- 
tion of  the  project  area  and  the  areas  immediately 
around  Bozeman  (Figure  20). 

Wetlands  with  the  potential  for  high  wetland 
hydrologic  functions,  including  water  storage, 
streamflow  maintenance,  and  groundwater  re- 
charge showed  variable  results  across  the  project 
area.  A  relatively  low  percentage  of  wetlands  had 
the  potential  for  high  functional  performance  of 
streamflow  maintenance,  with  the  majority  of  those 
wetlands  occurring  along  the  Gallatin  and  Jefferson 
Rivers  (Figure  21).  In  contrast,  the  majority  of  hy- 
drologic units  in  the  project  area  had  a  moderate  to 
high  percentage  of  wetlands  with  the  potential  for 
high  functional  performance  of  groundwater  re- 


charge (Figure  22).  The  west  slope  of  the  Bridger 
Mountains  had  the  highest  percentage  of  wetlands 
with  the  potential  for  high  water  storage  perfor- 
mance (Figure  23). 

Two  of  the  functions  related  to  biogeochemical 
processes,  shore  stabilization  and  nutrient  cycling, 
had  a  relatively  low  percentage  of  wetlands  across 
the  project  area  with  the  potential  for  high  func- 
tional performance  (Figure  24  and  Figure  25). 
Conversely,  many  of  the  hydrologic  units,  particu- 
larly those  along  the  East  Gallatin  and  JeflFerson 
Rivers  had  a  moderate  to  high  percentage  of  wet- 
lands with  the  potential  for  high  sediment  retention 
performance  (Figure  26). 

Functions  related  to  wildlife  and  plant  habitat, 
including  conservation  of  wetland  biodiversity, 
plant  community  maintenance,  terrestrial  habitat, 
and  aquatic  habitat,  displayed  interesting  patterns 
across  the  project  area.  The  river  valleys  had  the 
lowest  percentage  of  wetlands  with  the  potential 
for  high  conservation  of  wetland  biodiversity,  lack- 
ing ecologically  significant  wetland  types  such  as 
fens  and  forested  wetlands  (Figure  27).  In  contrast, 
these  areas  had  a  relatively  high  percentage  of  wet- 
lands with  the  potential  for  high  plant  community 
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Acres  of  Mapped  Wetlands 

Q"     0-179 
^1  180-501 
Bll  502-1.327 
^B  1.328-2.974 
^H  2.975-4.666 

Figure  1 7.    Wetland  landscape  profile  for  each  sixth  code  hydrologic  unit  across  the  Gallatin  project 
area:  wetland  acreage. 


Acres  of  Altered  Wetlands 

'  '  "    0-16 
■I  17- S7 
^B  58-153 
^H  154-373 


^H  374  -  579 

Figure  IS.   Wedand  landscape  profile  for  each  sixth  code  hydrologic  unit  across  the  Gallatin  project 
area:  acreage  of  altered  wetlands. 
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Acres  of  Isolated  Wetlands 

L 0-22 

IB  23-6-[ 
IB  62-110 
H  111  -  196 
^H  197-419 

Figure  19.    Wetland  landscape  profile  for  each  sixth  code  hydrologic  unit  across  the  Gallatin  project 
area:    acreage  of  isolated  wetlands. 


Percent  ol  Mapped  Wetlands  Located  on  Private  Land 
I  0-10 

^B  ''^-^ 

^■66-86 

^B  87-95 
■I  96-100 

Figure  20.   Wetland  landscape  profile  for  each  sixth  code  hydrologic  unit  across  the  Gallatin  project 
area:  acreage  of  wetlands  on  privately  owned  lands. 
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Percent  of  Wetlands  with  High  Streamflow  Maintenance  Functional  Pertormance 

I  00-0  1 

^H  0.2-0.4 
BB  0.5-0.9 
^B  1.0-8.2 
^B  8.3-16.2 

Figure  21.    H  elland  landscape  profile  of  wetlands  perfiarming  hydrologic  fitnclions  for  each  sixth  code 
hydrologic  unit  across  tlie  Gallatin  project  area:  streamfiow  maintenance. 


Percent  of  Wetlands  with  High  Groundwater  Recharge  Functional  Performance 

1 !  0-  16 

^H  17-35 
^B  36-56 
^B  57-79 
^B  80-100 

Figure  22.   Wetland  landscape  profile  of  wetlands  performing  hydrologic  functions  for  each  sixth  code 
hydrologic  unit  across  the  Gallatin  project  area:  groundwater  recharge. 
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Percent  of  Wetlands  with  High  Water  Storage  Functional  Performance 

■I  5-14 
■i  15-27 

^H  39-66 


Figure  23.   Wetland  landscape  profile  of  wetlands  performing  hydrologic  functions  for  each  sixth  code 
hydrologic  unit  across  the  Gallatin  project  area:  water  storage. 
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Percent  of  Wetlands  with  High  Shoreline  Stabilization  Functional  Performance 

I  I  0.00 

Ho.OI  -0  02 

^B  004-008 
^1  009-4  27 

Figure  24.   Wetland  landscape  profile  of  wetlands  performing  biogeochemical  functions  for  each  sixth 
code  hydrologic  unit  across  the  Gallatin  project  area:  shore  stabilization. 
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Percent  of  Wetlands  with  High  Nutrient  Cycling  Functional  Performance 

I  ^0-6 

^■7-13 
1^1  14-32 
IB  33-64 
^B  65-100 

Figure  25.    Wetland  landscape  profile  of  wetlands  performing  biogeochemical  functions  for  each  sixth 
code  hydrologic  unit  across  the  Gallatin  project  area    nutrient  cycling. 


Percent  of  Wetlands  with  High  Sediment  Retention  Functional  Performance 

[^0-7 
^H  8-22 
1^1  23    43 
Bi'>')-64 
■I  65    89 

Figure  26.   Wetland  landscape  profile  of  wetlands  performing  biogeochemical  functions  for  each  sixth 
code  hydrologic  unit  across  the  Gallatin  project  area:  sediment  retention. 
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maintenance,  particularly  along  the  Jefferson  River      habitat  function  were  concentrated  along  the  East 


(Figure  28).  Wetlands  in  these  areas  are  gener- 
ally wetter  and  thus  less  susceptible  to  invasion  by 
non-native  species  common  in  drier  wetland  types. 
Wetlands  with  the  potential  for  high  terrestrial 


Gallatin  River  (Figure  29),  whereas  most  wetlands 
with  the  potential  for  high  aquatic  habitat  function 
were  concentrated  in  areas  along  the  Madison  and 
portions  of  the  Jefferson  Rivers  (Figure  30). 


Percent  of  Wetlands  with  High  Conservation  of  Wetland  Biodiversity  Functional  Performance 

,  0-2 

H  3-6 
^H  7-13 


Figure  2  7.   Wetland  landscape  profile  of  wetlands  performing  wildlife  and  plant  habitat  functions  for 
each  sixth  code  hydrologic  unit  across  the  Gallatin  project  area:  conservation  of  biodiversity. 
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Percent  of  Wetlands  with  High  Native  Plant  Community  Maintenance  Functional  Performance 
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Figure  28.   Wetland  landscape  profile  of  wetlands  performing  wildlife  and  plant  habitat  functions  for 
each  sixth  code  hydrologic  unit  across  the  Gallatin  project  area:  plant  community  maintenance. 


Percent  of  Wetlands  with  High  Terrestrial  Habitat  Functional  Performance 
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Figure  29.   Wetland  landscape  profile  of  wetlands  performing  wildlife  and  plant  habitat  functions  for 
each  sixth  code  hydrologic  unit  across  the  Gallatin  project  area:  terrestrial  habitat. 
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Percent  of  Wetlands  with  High  Aquatic  Habitat  Functional  Perfonnance 
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Figure  30.   Wetland  landscape  profile  of  wetlands  performing  wildlife  and  plant  habitat  functions  for 
each  sixth  code  hydrologic  unit  across  the  Gallatin  project  area:  aquatic  habitat. 
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Discussion 


Our  analysis  shows  an  overall  increase  in  wetland 
area  based  on  mapping.  We  also  found  a  change 
in  wetland  types  Many  wetlands  mapped  as  palus- 
trine  scrub  shrub  in  the  1980"s  mapping  are  now 
palustrine  emergent  wetlands.  E.xamination  of  the 
aerial  imagery  revealed  that  much  of  the  change  is 
attributable  to  agricultural  modifications  (e.g..  live- 
stock grazing,  stream  dewatering.  and  conversion 
to  hay  pasture),  although  some  change  may  be  at- 
tributable to  the  mapping  itself  (e.g.  areas  mapped 
as  wetlands  in  the  1980's  being  mapped  as  riparian 
now). 

Despite  overall  wetland  increases  across  the  entire 
project  area,  concentrated  wetland  losses  occurred 
in  a  few  areas.  In  particular,  much  of  the  wetland 
change  in  the  areas  immediately  around  Bozeman 
was  attributable  to  urban  and  rural  development. 
This  area  has  seen  rapid  growth  and  much  of  the 
valley  bottom  along  the  East  Gallatin  River  has 
been  subdivided.  The  wetland  landscape  profiling 
also  revealed  that  the  areas  around  Bozeman  con- 
tain wetlands  with  the  potential  for  high  functional 
perfomiance  of  several  wetland  functions,  includ- 
ing groundwater  recharge,  streamflow  maintenance, 
water  storage,  sediment  retention,  and  terrestrial 
habitat.  Continued  impacts  to  wetlands  in  these 
important  areas  will  reduce  the  ability  of  wetlands 
to  perform  these  functions,  potentially  resulting  in 
ecological  and  economic  losses  in  these  areas. 

It  is  also  important  to  note  that  differences  between 
the  scales  of  the  imagery  used  in  the  historic  and 
new  mapping  products  make  accurate  quantifica- 
tion of  wetland  change  problematic.  We  mapped 
over  five  times  more  lotic  wetland  acres  using  the 
2005  imagery.  The  historic  wetland  mapping  was 
digitized  at  a  1 :58.000  scale  and  exhibits  consider- 
ably more  spatial  error  than  the  current  mapping 
digitized  at  a  1 :  12.000  scale.  The  increased  preci- 
sion made  possible  by  larger  scale  imagery  can 
result  in  an  increase  in  wetland  acres  mapped.  To 
accurately  compare  historic  and  new  wetland  map- 
ping products,  one  would  need  to  map  wetlands 
at  a  scale  of  1 :58,000  to  correspond  to  the  scale 
at  which  the  historic  mapping  was  completed,  but 
current  mapping  standards  and  statewide  mapping 


goals  require  a  1 :  12.000  scale.  Furthermore,  limi- 
tations exist  with  any  wetland  mapping  effort  de- 
rived almost  exclusively  from  photointerpretation 
techniques  (Tiner  2005).  Factors  such  as  photo 
quality,  scale,  and  environmental  conditions  at  the 
time  of  photo  acquisition  can  affect  mapping  accu- 
racy. Digital  wetland  maps  are  static  and  may  not 
reflect  the  dynamic  nature  of  wetlands  subject  to 
drastic  annual  and  seasonal  fluctuations  in  size  and 
distribution.  Comparisons  between  historic  and  up- 
dated wetland  mapping  are  also  limited  by  differ- 
ences in  imagery  and  mapping  scale.  Technologi- 
cal advances  associated  with  mapping  techniques 
have  greatly  increased  our  capacity  to  capture  more 
detailed  infonnation.  Additionally,  temporarily  and 
seasonally  flooded  wetlands  that  typically  dry  out 
by  the  end  of  the  growing  season  may  be  difficult 
to  reliably  separate  from  upland  areas.  We  also 
note  that  the  wetland  mapping  used  in  this  report 
was  completed  based  on  2005  imagery,  so  there  is  a 
nearly  five  year  time  lapse,  and  very  recent  changes 
are  not  captured. 

Limitations  also  exist  with  wetland  functional  as- 
sessments based  solely  on  photointerpretation  and 
best  professional  judgment.  We  emphasize  that 
the  functional  capacity  ratings  assigned  to  wet- 
lands in  this  project  are  only  potential  capacities. 
Data  on  actual  functional  capacity  would  require 
extensive  field  checking.  The  estimated  changes 
in  functional  capacities  may  also  be,  in  part,  a 
by-product  of  changes  in  mapping  conventions, 
specifically  the  addition  of  the  riparian  classifica- 
tion system.  Palustrine  forested  and  palustrine 
scrub-shrub  wetlands  associated  with  lotic  features, 
which  are  most  effective  at  performing  hydrologic 
and  habitat  functions,  may  have  been  overmapped 
in  the  1980"s  NWI  wetland  mapping.  Our  wetland 
change  analysis  revealed  that  some  wetlands  origi- 
nally mapped  as  palustrine  forested  and  palustrine 
scrub-shrub  wetlands  were  mapped  as  riparian  for- 
est and  riparian  scrub-shrub  in  the  new  mapping, 
so  this  reduction  in  mapped  wetland  acres  influ- 
ences our  estimated  changes  in  functional  capac- 
ity. Because  riparian  areas  were  not  mapped  in  the 
original  mapping  effort,  many  forested  and  scrub- 
shrub  areas  along  streams  and  rivers  were  mapped 
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as  wetland  rather  than  being  excluded  completely. 
Additionally,  it  is  often  difficult  to  distinguish  be- 
tween palustrine  and  riparian  types  based  on  aerial 
imagery  alone,  and  the  original  wetland  mapping 
had  only  cursory  ground-truthing.  In  addition  to 
differences  in  photointerpretation,  the  dynamic  na- 
ture of  lotic  systems  naturally  influences  the  type, 
location,  and  extent  of  ffoodplain  wetlands;  and 
changes  can  be  expected  over  multiple  decades. 
However,  anthropogenic  factors  such  as  dams,  rip- 
rap, and  increased  water  demands  have  almost  cer- 
tainly impacted  the  hydrology  of  these  floodplain 
wetlands. 

Additionally,  this  assessment  does  not  consider  the 
condition  of  the  adjacent  uplands  or  the  wetland 
itself;  two  factors  critical  to  assessing  the  integrity 
of  the  wetland.  Wetland  function  depends  largely 
on  wetland  type  (Brinson  1993),  but  we  could  not 
evaluate  the  effect  that  changes  in  wetland  type 
between  the  mid  to  late  1980"s  and  2005  have  had 
on  changes  in  wetland  function.  Thus,  extensive 
ecological  integrity  assessments  are  required  to 
fully  capture  the  condition  and  functional  capacity 
of  wetlands  in  the  study  area. 


This  analysis  should  be  considered  a  preliminary 
assessment  of  changes  in  the  Gallatin  Valley  wet- 
lands and  wetland  functional  capacity.  Although 
this  analysis  provides  general  information  on  the 
ability  of  wetlands  in  the  area  to  perform  certain 
functions,  it  does  not  assess  differences  in  func- 
tional capacity  among  wetlands  of  similar  type  and 
function.  Despite  these  limitations,  data  from  this 
analysis  can  provide  very  effective  conservation 
tools.  For  example,  detailed  wetland  and  riparian 
maps  can  be  used  to  identify  areas  with  the  poten- 
tial to  perform  wetland  functions  most  eflFectively, 
allowing  natural  resource  managers  and  other 
stakeholders  to  focus  or  prioritize  their  conserva- 
tion and  restoration  efforts.  The  wetland  landscape 
profiling  provides  an  understanding  of  the  role  of 
wetlands  across  the  project  area,  particularly  in  in- 
creasingly developed  landscapes.  Finally,  accurate 
digital  wetland  data  at  a  watershed  scale  can  also 
provide  the  foundation  for  initiation  and  tracking 
of  watershed  restoration  and  monitoring  efforts,  as 
well  as  a  tool  to  inform  watershed  councils,  local 
planners,  land  managers,  and  other  stakeholders 
about  the  types  and  location  of  wetlands  within  the 
Gallatin  Valley  and  the  fijnctions  they  perform. 
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Appendix  A.  Classification  of  Wetlands  and  Deepwater 
Habitats  of  the  United  States  (Cowardin  et  al.  1979) 
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Appendix  B.  Classification  of  Riparian  Systems 

(USFWS  1997) 
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Appendix  C.  Flowchart  for  Hydrogeomorphic  (HGM) 
Coding  of  Wetland  Polygons 
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Appendix  D.  Key  to  Hydrogeomorphic  (HGM)  Modifiers, 

Adapted  from  Tiner  (2003) 


Key  A-1:  Key  to  Wetland  Landscape  Position 

1.    Wetland  is  completely  surrounded  by  upland  (non-hydric  soils) Terrene 

Go  to  couplet  "d'  below 

1.  Wetland  is  not  surrounded  by  upland  but  is  connect  to  a  waterbody  or  other  wetland 2 

2.  Wetland  is  located  in  or  along  a  lake  or  reservoir  (permanent  waterbody  where 
standing  water  is  typically  deeper  than  6.6  feet  or  larger  than  20  acres),  including 
streamside  wetlands  in  a  lake  basin Lentic 

Go  to  couplet  'a'  below 
Note:  Lentic  wetlands  consist  of  all  wetlands  in  a  lake  basin  (i.e.,  the  depression  contain- 
ing the  lake),  including  lakeside  wetlands  intersected  by  streams  emptying  into  the  lake. 
The  upstream  limit  of  lentic  wetlands  is  defined  by  the  upstream  influence  of  the  lake, 
which  is  usually  approximated  by  the  limits  of  the  lake  basin.  Streamside  lentic  wetlands 
are  designated  Throughflow,  thereby  emphasizing  the  stream  flow  tlirough  these  wet- 
lands. 

2.  Wetland  does  not  occur  along  a  lake  or  reservoir 3 

3.  Wetland  is  located  in  a  river  or  stream  (including  in-stream  ponds),  within  its  banks, 

or  on  its  floodplain  and  is  periodically  flooded  by  the  river  or  stream 4 

3.  Wetland  is  not  located  in  a  river  or  stream  or  on  its  floodplain  OR  wetlands  along 

streams  are  NOT  subject  to  frequent  overflows Terrene 

4.  Wetland  is  the  source  of  a  river  or  stream,  but  there  is  no  channelized  flow  through 

the  wetland Terrene 

4.  Wetland  is  located  in  a  river  or  stream,  within  its  banks,  or  on  its  floodplain 5 

5.  Wetland  is  associated  with  a  river  (mapped  as  a  2-lined  watercourse  on  1 :24.000 

U.S.G.S.  topographic  map)  or  its  floodplain Lotic  River 

Go  to  Modifiers  "b".  'c\  and  'd'  below 
5.    Wetland  is  associated  with  a  stream  (single-line  watercourse  on  1 :24.000  U.S.G.S. 

topographic  map)  or  its  floodplain Lotic  Stream 

Go  to  Modifiers  'b'.  'c\  and  'd*  below 

a.    Wetland  is  associated  with  a  natural  water  body Natural  Lake 

a.    Wetland  is  associated  with  a  waterbody  created  by  a  dammed 

river  valley Dammed  River  Valley 

a.    Wetland  is  associated  with  a  waterbody  created  by  another 

obstruction Other  Dammed  Lake 

a.    Wetland  is  associated  with  a  waterbody  that  has  been  excavated Excavated  Lake 
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Slope  Modifiers 

b.    Water  flow  is  generally  slow  with  extensive  floodplain  development  (Cowardin's 

Lower  Perennial  Subsystem) Low  Gradient 

b.    Water  flow  is  moderate  to  fast  with  little  to  some  floodplain  development 

(Cowardin's  Upper  Perennial  Subsystem) Middle  Gradient 

b.  Water  flow  is  rapid  due  to  steep  gradient  with  little  or  no  floodplain 

development  (Cowardin's  Upper  Perennial  Subsystem) High  Gradient 

Water  Flow  Modifier 

c.  Water  flow  is  intermittent  or  ephemeral  (as  indicated  on  the  high  resolution 

National  Hydrography  Dataset) Intermittent 

Other  Special  Modifiers 

d.  Wetland  is  modified  by  a  beaver  dam Beaver 

d.    Wetland  is  modified  by  a  ditch  or  has  been  partially  drained Partially  Drained/Ditch 

d.    Wetland  has  been  farmed Farmed 

d.    Water  flow  or  wetland  modified  by  dam.  dike,  or 

impoundment Dammed/Diked/Impounded 

d.    Wetland  or  waterbody  has  been  modified  by  excavation Excavated 

d.    Wetland  forms  a  waterbody  less  than  6.6  feet  deep  or  smaller  than  20  acres. 

This  includes  in-stream  ponds Pond 

KeyB-1:  Key  to  Landforms 

I.    Wetland  hydrology  is  largely  influenced  by  groundwater  discharge  to  the  surface 

(can  occur  on  nearly  flat  landscapes)  or  wetland  occurs  on  a  slope  of  at  least  4% Slope 

1.  Wetland  hydrology  is  influenced  by  a  number  of  sources 2 

2.  Wetland  occurs  on  an  island Island 

2.  Wetland  does  not  occur  on  an  island 3 

3.  Wetland  occurs  with  the  banks  of  a  river  or  stream  or  along  the  shores  of  a  pond.  lake, 
or  island  and  is  either  ( 1 )  vegetated  and  at  least  saturated  or  semi-permanently  flooded 

or  (2)  a  non-vegetated  bank  or  shore  that  is  temporarily  or  seasonally  flooded Fringe 

3.  Wetland  does  not  occur  along  a  bank  or  shore 4 

4.  Wetland  occurs  on  a  floodplain  of  a  Lotic  River Floodplain 

4.  Wetland  does  not  occur  on  a  floodplain  of  a  Lotic  River 5 

5.  Wetland  exists  in  a  distinct  depression  in  various  landscape  positions Basin 

5.    Wetland  receives  virtually  no  groundwater  discharge  (playa) Flat 
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Key  C-1:  Key  to  Water  Flow  Paths 

1.    Water  levels  fluctuate  due  to  lake  influences  or  variable  river  levels,  but  water  does 

not  flow  through  the  wetland Bidirectional 

1.  Wetland  is  not  influenced  by  fluctuating  water  levels 2 

2.  Wetland  receives  surface  or  groundwater  from  a  watercourse,  other  waterbody,  or 
from  another  wetland  at  a  higher  elevation  and  surface  or  groundwater  passes  through 

it  to  a  watercourse,  other  waterbody,  or  another  wetland  at  a  lower  elevation Throughflow 

2.  Water  does  not  pass  through  the  wetland  to  other  wetlands  or  waters 3 

3.  Wetland  is  in  a  closed  depression  (based  on  closed  elevation  contours)  or  flat  areas 
where  water  conies  from  surface  water  runoff  and/or  groundwater  discharge.  Isolated 
wetlands  lack  channelized  surface  water  inflow  and  outflow Isolated 

3.  Wetland  is  not  hydrologically  or  geographically  isolated 4 

4.  Wetland  exists  in  a  sink  where  no  outlet  exists,  yet  water  enters  via  an  intermittent 

or  perennial  stream  or  from  a  wetland  at  a  higher  elevation Inflow 

4.    Wetland  receives  no  surface  or  groundwater  inflow  from  a  wetland  or  other 
permanent  waterbody  at  a  higher  elevation  and  water  is  discharged  from  the 
wetland  to  a  watercourse,  waterbody,  or  other  wetland  at  a  lower  elevation Outflow 

KeyD-1:  Key  to  Waterbody  Types 

1.    Waterbody  is  predominantly  flowing  water 2 

1.  Waterbody  is  predominantly  standing  water 3 

2.  Waterbody  is  a  polygonal  feature  on  a  U.S.G.S.  topographic  map  (1 :24.000) River 

2.  Waterbody  is  a  linear  feature  on  a  U.S.G.S.  topographic  map  ( 1 :24.000) Stream 

Go  to  couplets  "b\  'c',  and  'd'  below 

3.  Waterbody  is  permanently  flooded  and  deep  (>  6.6  feet) Deep  Water 

Go  to  couplet  'a'  below 

a.  Waterbody  is  created  by  a  dammed  river  valley Dammed  River  Valley 

a.  Waterbody  is  created  by  another  obstruction Other  Dammed  Lake 

a.  Waterbody  has  been  excavated Excavated  Lake 

a.  Waterbody  is  natural Natural  Lake 
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Slope  Modifiers 

b.    Water  flow  is  generally  slow  with  extensive  floodplain  development.  Stream 

slope  is  less  than  2%.  (Cowardin's  Lower  Perennial  Subsystem) Low  Gradient 

b.    Water  flow  is  moderate  to  fast  with  little  to  some  floodplain  development. 
Stream  slope  is  between  2  and  4%.  (Cowardin's  Upper  Perennial 
Subsystem) Middle  Gradient 

b.  Water  flow  is  rapid  due  to  steep  gradient  with  little  or  no  floodplain  development. 
Stream  slope  is  greater  than  4%.  (Cowardin's  Upper  Perennial 

Subsystem) High  Gradient 

Water  Flow  Modifier 

c.  Water  flow  is  intermittent  or  ephemeral  (as  indicated  on  the  high  resolution 

National  Hydrography  Dataset) Intermittent 

Other  Special  Modifiers 

d.  Water  flow  is  dammed Dammed 
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Appendix  E.  Procedure  for  Coding  Wetland  Polygons  with 
Hydrogeomorphic  Modifiers  in  ArcGIS 


The  hydrogeomorphic  (HGM)  coding  process  has  been  semi-automated  to  populate  each  wetland 
polygon  with  an  HGM  code.  Preparation  for  the  semi-automated  process  involves  exporting  the 
wetland  polygons  to  a  new  layer  and  adding  fields  for  each  of  the  HGM  code  components.  Once 
the  HGM  coding  process  has  been  completed  the  exported  table  can  be  joined  to  the  original 
NWI  geodatabase  to  populate  the  geodatabase  HGMCODE  field. 

1 .  The  average  slope  is  calculated  for  each  of  the  wetland  polygons.  Then  the  following  fields 
are  added  to  the  table:  Waterbody,  Gradient.  LakeMod.  WaterFlow.  Pond,  SpecMod, 
Landform.  and  FlowPath. 

2.  Each  field  is  then  populated  in  succession.  Seven  waterbody/landscape  position  types  are 
included  in  the  HGM  code:  deep  water  (DW).  lentic  (LE),  riverine  (RV),  stream  (ST),  lotic 
river  (LR),  lotic  stream  (LS),  and  terrene  (TE).  The  "DW"  code  is  assigned  to  all  wetlands 
attributed  as  an  LI .  Those  wetlands  coded  as  L2  are  given  the  "LE"  code.  Additionally,  all 
wetlands  within  40  meters  of  lakes  with  an  average  slope  less  then  4%  are  coded  as  "LE". 
The  "RV"  code  is  assigned  to  rivers  attributed  as  R2UBH  or  R3UBH.  All  other  streams  are 
coded  as  "ST"  including  R2/3UBF,  R2/3UBG,  and  R4SB  streams.  Wetlands  coded  as  "RV" 
within  300  meters  of  rivers  and  with  an  average  slope  less  than  4%  are  assigned  the  "LR" 
code.  Wetlands  coded  as  "ST"  within  100  meters  of  rivers  and  an  average  slope  of  less  that 
4%  are  coded  as  "LS".  Then  using  the  high  resolution  national  hydrography  dataset  (NHD) 
layer,  the  perennial  and  intermittent  streams  are  buffered  by  1 00  meters  and  the  ephemeral 
streams  are  buffered  by  20  meters.  All  wetlands  that  intersect  either  the  100  meter  or  20  me- 
ter buffer  with  an  average  slope  of  less  than  4%  are  also  assigned  the  "LS"  code.  The  "TE" 
code  is  then  assigned  to  any  remaining  uncoded  wetlands. 

3.  Using  the  average  slope  value,  each  river  and  stream  is  assigned  a  gradient  value.  Those 
streams  with  an  average  slope  less  than  2%  are  given  a  "1"  or  low  gradient  value.  Streams 
with  an  average  slope  between  2%  and  4%  are  assigned  a  value  of  "2".  A  value  of  "3"  is  as- 
signed those  streams  with  an  average  slope  greater  than  4%.  Additionally,  streams  attributed 
as  R4SB  or  with  a  water  regime  of  F  (semipermantently  flooded)  are  assigned  a  value  of  "4" 
for  intermittent  in  the  WaterFlow  field. 

4.  Each  deep  water  and  lentic  wetland  is  assigned  a  lake  modifier.  Natural  lakes  are  given  a 
"1".  Dammed  lakes  are  assigned  a  "3".  And  excavated  lakes  are  assigned  the  value  "4". 
The  lentic  wetlands  are  given  the  same  lake  modifier  as  the  lake  with  which  they  are  associ- 
ated. 

5.  The  Spec_mod  field  is  populated  with  the  special  modifiers  that  a  wetland  polygon  may  con- 
tain. The  possible  special  modifiers  are  "b"  (beaver),  "h"  (dammed/impounded),  "x"  (exca- 
vated), "d"  (drained),  and  "f '  (farmed).  Not  all  wetlands  will  have  a  special  modifier. 

6.  The  Pond  field  is  populated  with  a  "p"  for  all  wetland  polygons  attributed  as  PAB  (palustrine 
aquatic  bed)  or  PUB  (palustrine  unconsolidated  bottom). 
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7.  The  Landform  field  is  populated  for  each  waterbody  type.  There  are  three  possible  landform 
values  possible  for  the  lentic  system.  Wetlands  that  are  islands  within  lakes  are  assigned 
"IL".  Wetlands  with  the  water  regimes  A  or  C  and  ponds  are  "BA"  or  basins.  Wetlands 
with  the  water  regimes  B  or  F  and  shores  along  rivers  and  lakes  are  assigned  "FR"  or  fringe. 
Ponds  near  lotic  rivers  and  streams  are  also  basins  ("BA").  All  other  wetlands  along  lotic 
streams  are  also  assigned  "BA"  except  those  wetlands  with  the  water  regimes  B  or  F  and 
shores,  which  are  defined  as  "FR".  Along  lotic  rivers,  those  wetlands  having  water  regimes 
B  or  F  and  shores  are  also  defined  as  "FR".  but  all  other  wetlands  are  defined  as  "FP"  or 
floodplain.  All  terrene  wetlands  are  defined  as  basins  ("BA")  except  those  with  an  average 
slope  greater  than  4%,  which  are  defined  as  slope  or  "SL". 

8.  The  Flowpath  field  is  populated  for  each  wetland  type.  All  lotic  rivers  and  streams  are 
populated  with  "TH"  or  throughflow.  Wetlands  attributed  with  L2  a  water  regime  of  B  or  F 
or  shores,  (excluding  ponds)  are  attributed  as  "BI"  or  bidirectional.  All  other  lentic  wetlands 
are  defined  as  throughflow  or  "TH".  All  dammed  lakes  are  also  defined  as  "TH".  Lakes  or 
"DW"  can  also  be  defined  as  "IN".  "OU",  or  "IS"  for  inflow,  outflow,  and  isolated.  Each 
lake  is  checked  against  the  digital  USGS  topographic  map  to  determine  the  Flowpath.  Ter- 
rene wetlands  are  defined  as  either  "IS"  (isolated)  or  "CO"  (complex).  All  terrene  wetlands 
that  are  greater  than  20  meters  from  another  wetland  are  defined  as  isolated. 

9.  Once  each  field  has  been  populated,  the  HGM_CODE  field  can  be  updated  using  the  field 
calculator.  The  deep  water  code  consists  of  the  Waterbody,  LakeMod.  and  a  Flow_Path. 
The  riverine  code  contains  the  Waterbody,  Gradient,  and  SpecMod  fields.  The  stream  code 
includes  the  fields  Waterbody.  Gradient,  Water_Flow,  and  Spec_Mod.  Lentic  polygons 
include  the  Waterbody.  Lake_Mod,  Landform,  and  Flow_path  fields.  Finally  the  lotic  rivers 
and  streams  and  terrene  wetlands  contain  the  fields:  Waterbody,  Pond,  SpecMod.  Land- 
form,  and  Flow_path. 
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Appendix  F.  Relative  Functional  Performance  Levels  for 

Wetlands  Classified  with  National  Wetland  Inventory 

(NWI)  AND  Hydrogeomorphic  (HGM)  Codes. 

(A  RANK  OF  "1"  INDICATES  A  HIGH  PERFORMANCE  LEVEL; 

A  RANK  OF  "2"  INDICATES  A  MODERATE  PERFORMANCE  LEVEL; 

AND  A  RANK  OF  "3"  INDICATES  A  LOW  PERFORMANCE  LEVEL) 
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2 

PEMAd 

LSdBATH 

22  4 

2 

2 

2 

2 

2 

PEMAI' 

LRfFPTH 

5'I9  1 

2 

2 

2 

2 

2 

PEMA/ 

LSfBATH 

517  7 

2 

2 

2 

2 

2 

PEMAf 

TEIBACO 

522 

1 

2 

2 

2 

PEMAf 

TEffiAIS 

43  2 

1 

2 

2 

2 

PEMAf 

TEfSLCO 

08 

1 

2 

2 

2 

PEMAI" 

TEISLIS 

14  11 

1 

2 

2 

2 

PEMAh 

le:batfi 

55  5 

1 

3 

1 

2 

PEMAh 

TEhBAIS 

0  1 

1 

3 

2 

2 

PEMAh 

TEhSLCO 

6  1 

1 

3 

2 

2 

PEM,\h 

TEbSLIS 

35 

1 

3 

2 

2 

PEMAh 

LRbFPTH 

1  4 

t 

1 

2 

2 

PEMAh 

LShBATH 

('2  1 

1 

1 

1 

2 

PEMA\ 

TExBACO 

128 

1 

3 

2 

2 

PEMAn 

TExBAIS 

19  8 

1 

3 

2 

2 

PEMAx 

TExSLCO 

43 

1 

3 

2 

2 

PEMA\ 

TExSLIS 

38 

1 

3 

2 

PEMAx 

LRxFPTH 

35  5 

1 

1 

2 

PEMAx 

LSxBATH 

1906 

1 

1 

2 

PEMB 

TESLIS 

58 

1 

2 

2 

PFMB 

LSFRTH 

14  6 

2 

1 

3 

PtMC 

LRFPTH 

2609  2 

2 

3 

3 

PEMC 

TEBAIS 

23  1 

1 

3 

2 

PEMC 

TESLIS 

30  8 

1 

3 

2 

PEMC 

TEBACO 

75  3 

I 

3 

2 

3 

PEMC 

TESLCO 

27  3 

1 

3 

2 

3 

PEMC 

LSBATH 

2598  6 

2 

2 

2 

2 

PEMCd 

LRdFPTH 

15  3 

2 

2 

2 

2 

PEMCf 

TEIBACO 

54 

1 

3 

2 

1 

PEMCr 

TEfBAIS 

3  3 

1 

3 

2 

' 

PEMCf 

LR/FPTH 

75  8 

2 

2 

2 

PEMCr 

LSfBATH 

33 

2 

2 

2 

PEMCh 

LRhFPTH 

1  5 

2 

3 

3 

1 

PEMCh 

LShBATH 

30  8 

2 

2 

3 

PEMCh 

LE2BATH 

33  5 

1 

1 

2 

PEMCh 

TEhBACO 

0  7 

I 

1 

1 

PEMCh 

TEhSLCO 

1  2 

1 

1 

1 

1 

PEMCh 

TEhSLIS 

1  1 

1 

1 

J 

1 
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NWI 
Code 

HGM 
Code 

Acres 

Water 
Storage 

SlreamlloH 
Maintenance 

(Groundwater 
Recharge 

Nutrient 
Cycling 

Sediment 
Retention 

Shoreline 
Stabilization 

Native  Plant 
Community 
Maintenance 

Terrestrial 
Habitat 

Aquatic 
Habitat 

Conserxalion 
ofWetland 
Biodiversity 

PEMC\ 

LRxFPTH 

75  0 

: 

3 

PEMC\ 

TExBACO 

28  3 

2 

2 

PEMCx 

TExBAIS 

22  2 

2 

2 

PEMCn 

TExSLCO 

26 

2 

2 

PEMCx 

TExSLIS 

III  1 

2 

2 

PEMC\ 

LSxBATH 

453  4 

1 

2 

PENtF 

TEBACO 

48 

3 

2 

PEMF 

TEBAIS 

1  ^ 

3 

2 

PEMF 

TESLIS 

1  I 

3 

2 

PEMF 

LRFRTH 

\M  1 

3 

3 

PEMF 

LSFRTH 

3f.:  (1 

3 

1 

PEMFb 

LRbFRTH 

0  5 

3 

3 

PEMFh 

LRhFRTFI 

1  5 

2 

3 

2 

PEMFh 

LShFRTH 

11)3 

2 

3 

2 

PEMFh 

TEhBACO 

0  5 

2 

1 

2 

PEMFh 

TEhBAIS 

0  3 

2 

1 

2 

PEMFh 

TEhSLlS 

0  2 

2 

1 

2 

PEMFx 

LRxFRFH 

C>4 

2 

3 

2 

PEMFx 

LSxFR™ 

'>0 

2 

3 

2 

PEMFv 

TExBACO 

06 

3 

2 

2 

PEMFn 

TExBAIS 

0  7 

3 

2 

2 

PEMK-X 

TExBACO 

48  8 

1 

1 

I 

PFOA 

l.SBATH 

3  } 

2 

3 

2 

2 

PFOA 

LRFPTH 

1  3 

2 

3 

2 

2 

PSSA 

LRFPTH 

23(.l  3 

2 

3 

2 

2 

2 

PSSA 

le:bath 

0  7 

2 

3 

2 

2 

2 

PSSA 

TEBACO 

127  1 

3 

3 

2 

2 

2 

PSSA 

TEBAIS 

52  5 

3 

2 

2 

2 

PSSA 

TESLIS 

1 70  4 

3 

2 

2 

1 

PSSA 

TESLCO 

378  6 

2 

3 

2 

1 

PSSA 

LSBATH 

21X)4  (i 

2 

2 

2 

2 

PSSAb 

LSbBATH 

0  5 

2 

2 

2 

2 

PSS,Ah 

TEhBACO 

0  2 

3 

2 

2 

2 

PSSAh 

TEhSLCO 

1  9 

3 

2 

2 

2 

PSSAh 

LRhFPTH 

04 

2 

1 

3 

1 

1 

PSS.\h 

LSbBATH 

3  7 

2 

1 

3 

1 

1 

PSSAx 

LRxFPTH 

5  7 

2 

2 

3 

2 

2 

PSSAx 

TExBACO 

44 

3 

2 

2 

2 

2 

PSSAx 

TExBAIS 

1  0 

3 

2 

2 

2 

2 

PSSAx 

TExSLCO 

1  3 

3 

2 

2 

2 

2 

PSSA^ 

TExSLIS 

0  6 

3 

2 

2 

2 

2 

PSSAx 

LSxBATH 

37  7 

1 

3 

1 

2 

1 

PSSC 

LRFPTH 

705  ^ 

2 

3 

2 

2 

PSSC 

TEBACO 

3  1 

3 

2 

2 

1 

PSSC 

TESLCO 

1  1 

3 

2 

2 

2 

PSSC 

TESLIS 

2  7 

3 

2 

2 

1 

PSSC 

LSBATH 

217  7 

2 

2 

2 

2 

PSSCb 

LSbBATH 

0  4 

2 

3 

2 

2 

PSSCh 

TEhBACO 

02 

3 

2 

2 

2 

PSSCx 

LRxFPTH 

3  3 

2 

3 

2 

2 

2 

PSSCx 

LSxBATH 

16  8 

3 

2 

3 

2 

2 

2 

PSSCx 

TExBACO 

0  6 

1 

2 

2 

1 

1 

PSSCx 

TExSLCO 

0  4 

3 

2 

2 

2 

PSSF 

LRFRTH 

2  6 

3 

3 

3 

3 

PSSF 

LSFRTH 

1  3 

3 

3 

3 

3 

PSSF 

TEBACO 

1  0 

3 

2 

2 

2 

PSSFb 

LRbFRTH 

27 

2 

3 

2 

2 

PSSFb 

LSbFRTH 

0  2 

2 

3 

2 

2 

PSSFb 

LRhFRTH 

0  1 

1 

3 

1 

PUBFh 

TEpbBAIS 

0  1 

1 

1 

1 

PUBFx 

LRpxBATH 

24  •» 

1 

2 

1 

PUBFx 

LSpxBATH 

35  4 

1 

2 

1 

PUBFx 

TEpxBACO 

59 

2 

2 

1 

PLIBFx 

TEpxBAIS 

183 

1 

1 

1 

PUBGx 

LSpxBATH 

94 

1 

1 

PUBGx 

TEpxBACO 

24  4 

2 

1 

PUBGn 

TEpxBAIS 

54 

1 

1 

PLIBRx 

TEpxBAIS 

0  7 

1 

1 

PUSA 

LSBATH 

2  9 

1 

2 

2 

PUSA 

LRFPTH 

2  5 

2 

2 

2 

PUSA 

TEBACO 

0  1 

2 

2 

1 

PUSA 

TEBAIS 

2  5 

2 

2 

1 

PUSAh 

LShBATH 

85 

1 

2 

1 

PUSAh 

le:frth 

3  6 

2 

2 

2 

PUS.Ah 

TEhBACO 

1  5 

2 

2 

1 

PUSAh 

TEbBAlS 

1  0 

2 

2 

1 

PUSA.X 

LSxBATH 

1  r. 

1 

2 

2 

2 
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NWI 

Code 

HGM 

Code 

Acres 

Water 
Storage 

Stream  flow 
Maintenance 

Groundwater 
Recharge 

Nutrient 
Cycling 

Sediment 
Retention 

Shoreline 
Stabilization 

Native  Plant 
Community 
Maintenance 

Terrestrial 
Habitat 

Aquatic 
Habitat 

Conservation 
orWetland 
Biodiversity 

PLSAx 

I  RxIPlH 

(14 

2 

2 

2 

2 

2 

PUSAv 

THvBACO 

Oil 

2 

1 

2 

PUSAx 

Tt\BAIS 

2  5 

2 

1 

2 

PUSC 

LSBAl-H 

5') 

1 

2 

2 

PUSC 

LRFPIH 

1114 

2 

2 

2 

PlSC 

TEBACO 

3  1 

2 

1 

2 

PUSC 

TtBAlS 

1 ; 

2 

1 

2 

PLiSCh 

LShBATH 

2  3 

1 

2 

2 

PUSCh 

LRhFPTH 

(14 

: 

2 

2 

PLSCh 

le;frth 

42 

2 

2 

2 

PUSCh 

TFhBACO 

(Ul 

2 

2 

PUSCh 

TthBAIS 

1  1 

2 

2 

PlSCx 

ISsBAFH 

1.  1 

1 

2 

PUSCx 

LRnFP™ 

2  5 

2 

2 

PUSC\ 

TE\BACO 

4  9 

2 

2 

PL'SC\ 

TExBAlS 

115 

1 

R2ABG 

STl 

22  4 

3 

R2ABG 

ST14 

12') 

3 

r;abg 

ST2 

72 

3 

R2.\BG 

ST24 

20 

3 

R2ABG 

ST3 

II  0 

3 

RJABG 

ST34 

0  1 

2 

3 

R2UBF 

STl 

>>3  •) 

2 

3 

R2UBF 

STI4 

513 

2 

3 

r;ubf 

ST2 

5  5 

2 

3 

R21IBH 

st:4 

76 

2 

3 

r:;i!BF\ 

STUk 

154 

2 

2 

R2UBFS 

sriN 

57  5 

2 

2 

R2UBFX 

ST2\ 

6  8 

2 

2 

R2llBF\ 

ST34X 

20  3 

2 

2 

r:ubo 

STl 

(.4  5 

2 

3 

R^l'BG 

STI4 

l.S(. 

2 

3 

R2UBG 

ST2 

1  5 

2 

3 

r;ubg 

ST3 

0  4 

2 

3 

R2UBH 

RVI 

4103  8 

2 

3 

R2USA 

LRFRTH 

083  2 

2 

3 

2 

2 

R2USA 

LSFRTH 

34 

2 

3 

2 

2 

r:;usc 

LSFRTH 

14 

2 

3 

2 

2 

R21ISC 

LRFRTH 

375  4 

3 

2 

RJL'BF 

STI4 

10  0 

2 

3 

RniBFh 

STMh 

0  6 

2 

2 

R3UBFN 

STI4\ 

6  1 

2 

2 

R-UIBF\ 

ST24S 

ISO 

2 

2 

R3UBF\ 

ST34\ 

54  4 

2 

2 

R3UBG 

STl 

0  0 

2 

3 

R3UBG 

ST14 

21  7 

2 

3 

R3LIBH 

RV3 

75  2 

2 

3 

R3UBH 

RVI 

182  5 

1 

3 

R3USA 

LRFRTH 

15  1 

2 

2 

R3USC 

LRFRTH 

1  8 

3 

2 

R4SBA 

STI4 

7  7 

1 

2 

R4SBA 

ST24 

79 

1 

2 

R-tSBA 

ST34 

2  1 

1 

2 

R4SBA\ 

STUx 

22  1 

1 

2 

R4SBAX 

ST24V 

21  8 

1 

2 

R4SBA\ 

ST34\ 

68  1 

1 

2 

R4SBC 

STI4 

448 

1 

2 

R4SBC 

ST24 

II  6 

1 

2 

R4SBC 

ST34 

166 

1 

2 

R4SBCN 

STI4X 

1109 

2 

1 

2 

R4SBC\ 

ST24S 

51  5 

2 

1 

1 

R4SBC\ 

ST34\ 

59  9 

2 

1 

2 
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Appendix  G.  Acres  of  major  land  cover  and  land  use  types 

WITHIN  EACH  SIXTH  CODE  HYDROLOGIC  UNIT  IN  THE  GaLLATIN 
PROJECT  AREA  FROM  THE  2009  MONTANA  LaND  CoVER  THEME 


Sixth  Code  Hydrologic  Unit 

Open 
Water 

Developed 
Area 

Pasture/Cultivated 
Cropland 

Native 
Vegetation 

Harvested 
Forest 

Antelope  Creek 

38.3 

584.5 

3,390.8 

26,832.3 

0.0 

Bear  Creek 

0.0 

142.3 

873.8 

11,858.3 

113.2 

Big  Bear  Creek 

2.4 

165.0 

893.1 

11.552,7 

1,220.7 

Boulder  River-Conrovv  Creek 

2.0 

681.6 

1,325.2 

33.832.1 

0.0 

Bozeman  Creek 

11.3 

3.143.5 

3.837.9 

25.721.4 

356.1 

Burnt  Creek 

2.2 

98.1 

121.0 

11,827.6 

0.0 

Cottonwood  Creek 

0.4 

106.5 

950.1 

17,578.7 

294.7 

Dry  Creek 

0.0 

0.0 

1.3 

4.8 

0.0 

Dr>'  Hollow  Creek 

22.2 

421.2 

3.488.7 

4,867.1 

0.0 

East  Gallatin  River-Belgrade 

20.9 

4,038.7 

12.098.2 

6,804.1 

0.0 

East  Gallatin  River-Bozeman 

23.4 

5.035.4 

7.539.8 

7,155.5 

0.0 

East  Gallatin  River-Kelly  Creek 

0.0 

2.557.1 

3.826.5 

10,003.3 

0.0 

East  Gallatin  River-Nixon  Gulch 

9.6 

873.1 

8.519.5 

19.522.9 

0.0 

Elk  Creek 

0.0 

460.1 

3.343.5 

24.055.9 

0.0 

Ennis  Lake 

0.0 

0.0 

0.0 

11.3 

0.0 

Fish  Creek 

6.4 

480.1 

2.685.4 

5.642.8 

0.0 

Gallatin  River-Central  Park 

58.0 

3,227.6 

26.917.9 

11.682.8 

0.0 

Gallatin  River-Gal  latin  Gateway 

171.7 

1.879.7 

17.620.5 

18.565.2 

0.2 

Gallatin  River-Manhattan 

20.2 

2.340.7 

11.811.1 

15.127.2 

0.0 

Gallatin  River-Missouri  Headwaters 

35.4 

719.9 

2.561.8 

16.502.5 

0.0 

Gallatin  River-Spanish  Creek 

68.5 

108.3 

0.0 

9.818.7 

176.8 

Gallatin  River- Wilson  Creek 

95.4 

257.1 

2,451.0 

15,289.6 

584.7 

Jackson  Creek 

0.0 

4.2 

541.8 

10,732.7 

896.7 

Jefferson  River-Cardwell 

154.3 

865.3 

6,452.3 

26,039.4 

0.0 

Jefferson  River-Mill  Creek 

54.0 

242.0 

1.039.9 

4,831.5 

0.0 

Jefferson  River-Sappington 

280.7 

1.033.2 

1.130.7 

28.308.1 

38.5 

Jefferson  River- Three  Forks 

336.5 

2.654.9 

9.159.1 

24.986.6 

0.0 

Jefferson  Slough 

89.6 

2.970.3 

2.487.0 

9.295.8 

0.0 

Little  Pipestone  Creek 

0.0 

524.0 

672.5 

4.700.1 

0.0 

Lower  Big  Pipestone  Creek 

7.8 

1,529.2 

3.647.0 

25.386.9 

60.5 

Lower  Bridger  Creek 

0.0 

71.6 

1,142.7 

11,532.2 

22.2 

Lower  Camp  Creek 

5.8 

1.211.6 

18,692.0 

13,926.1 

0.0 

Lower  Cherry  Creek 

0.0 

6.4 

149.0 

17,001.6 

0.0 

Lower  Dry  Creek 

0.9 

482.8 

5,218.9 

13,472.6 

0.0 

Lower  Hot  Springs  Creek 

3.6 

517.1 

804.6 

15,113.4 

0.0 

Lower  Hyalite  Creek 

25.8 

5,553.0 

27,485.0 

5,248.7 

0.0 

Lower  North  Meadow  Creek 

0.0 

0.0 

0.0 

1.053.5 

0.0 

Lower  South  Boulder  Creek 

0.2 

314.9 

2.330.5 

8.313.3 

0.0 

Lower  Whitetail  Creek 

0.0 

1.200.3 

2.020.2 

13,232.7 

0.0 

Madison  River-Bear  Trap  Canyon 

139.2 

0.0 

0.9 

14,062.0 

0.0 

Madison  River-Red  Mountain 

247.5 

406.8 

205.5 

12,859.9 

0.0 

Madison  River-Three  Forks 

342.3 

2,583.1 

7,717.5 

22,337.7 

0.0 

Madison  River-Willow  Springs 

174.1 

796.6 

3,900.8 

23,978.3 

0.0 

Middle  Cherrv  Creek 

0.0 

0.0 

0.0 

6,435.6 

0.0 
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Sixth  Code  Hydrologic  Unit 

Open 
Water 

Developed 
Area 

Pasture/Cultivated 
Cropland 

Native 
Vegetation 

Harvested 
Forest 

Middle  South  Boulder  Creek 

0.0 

0.0 

18.2 

12.502.5 

0.0 

Milligan  Creek 

2.0 

1.051.5 

757.3 

28,408.2 

0.0 

Missouri  River-Clarkston 

283.1 

9.1 

2.253.3 

32,709.5 

34.0 

Missouri  River- Trident 

225.1 

91.8 

312.9 

13,975.2 

0.0 

Mud  Spring  Gulch 

3.8 

1,392.4 

3.058.4 

30,878.3 

0.0 

North  Fork  Spanish  Creek 

0.0 

0.0 

188.1 

3.893.7 

0.0 

North  Fork  Willow  Creek 

34.9 

427.7 

4.017.8 

15.390.8 

2.2 

Norwegian  Creek 

423.4 

733.5 

1,336.1 

21.297.1 

0.0 

Pass  Creek 

0.0 

205.5 

2,230.2 

10.631.3 

91.8 

Piedmont  Swamp 

10.2 

1.051.5 

3,653.0 

11,438.2 

0.0 

Pole  Creek 

0.0 

0.0 

6.7 

11,900.1 

0.0 

Rattlesnake  Creek 

19.8 

1,400.6 

3,758.7 

31,356.7 

0.0 

Rey  Creek 

0.0 

435.9 

5,426.2 

11.627.4 

0.0 

Rocky  Creek 

0.0 

505.7 

526.9 

19,954.3 

747.7 

Ross  Creek 

3.} 

492.6 

5,997.3 

6.667.6 

0.0 

Smith  Creek 

12.5 

622.9 

7,762.0 

19,999.0 

0.0 

South  Fork  Willow  Creek 

0.9 

338.3 

2,960.3 

3.894.3 

0.0 

Spanish  Creek 

21.1 

0.2 

0.0 

8.967.6 

0.0 

Upper  Bridger  Creek 

0.0 

0.0 

575.8 

7.245.6 

586.0 

Upper  Camp  Creek 

0.0 

175.9 

3,221.4 

10.355.8 

0.0 

Upper  Dry  Creek 

0.0 

762.1 

4,185.0 

29.747.0 

0.4 

Upper  Hot  Springs  Creek 

0.0 

4.7 

43.8 

2.839.5 

0.0 

Upper  South  Boulder  Creek 

0.0 

0.0 

0.0 

1.514.1 

46.7 

Upper  South  Fork  Sixteenmile  Creek 

0.0 

0.0 

0.0 

3.476.0 

359.6 

Wann  Springs  Creek 

0.0 

1,940.6 

5.737.1 

20.436.4 

0.0 

Willow  Creek 

51.3 

329.5 

1,206.3 

10.633.1 

125.4 
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